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Abstract 


This  research  developed  a  methodology  for  assessment  of  the  exposure  of  a 
mature  com  crop  to  trichloroethylene  from  a  contaminated  aquifer.  The 
methodology  was  then  applied  to  the  case  of  Hill  AFB  to  determine  the  ability  of  the 
methodology  to  provide  information  about  a  specific  exposure.  There  are  currently 
several  locations  where  ground  water,  contaminated  with  TCE  originating  from  an 
Air  Force  facility,  has  the  potential  to  expose  food  crops.  Current  procedure  for 
estimating  the  risk  to  human  health  from  this  type  of  food  crop  exposure  is  focused 
on  sampling  the  food  product  and  lacks  the  information  to  predict  where  exposure 
problems  are  more  likely  to  exist.  A  review  of  the  potential  exposure  pathways  from 
the  aquifer  to  the  crop  was  conducted.  On  the  basis  of  this  review,  the  exposures  due 
to  soil  gas  and  irrigation  were  modeled.  Empirical  estimates  were  used  to 
£q)proximate  the  ejqiected  flux  of  soil  gas  vaporizing  directly  from  the  aquifer.  On 
the  basis  of  this  approximation,  the  exposure  in  the  air  of  the  crop  canopy  was 
mathematically  estimated.  Analytical  models  were  developed  to  simulate  the 
amount  of  the  contaminant  reaching  the  crop  from  two  different  means  of  irrigation. 
The  subsequent  exposure  once  the  contaminated  irrigation  water  had  reached  the 
crop  was  modeled  both  in  the  air  of  the  crop  canopy  and  the  soil  phase  near  the  root 
system.  The  methodology  provided  insights  into  ^\ilich  exposure  pathways  are  more 
important  than  others  and  revealed  some  to  be  negligible.  Information  was  also 
gained  regarding  the  environmental  parameters  that  are  most  influential  in 
determining  the  amount  of  exposure. 
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EXPOSURE  OF  A  FOOD  CROP  TO  TRICHLOROETHYL  ENE  FROM  A 

CONTAMINATED  AQUIFER 

I.  Introduction 

Trichloroethy!'  le  is  a  solvent  common  to  many  commercial  operations.  Its 
wide  use  has  included  everything  from  degreasing  fabricated  metal  parts  to  use  as  a 
general  anesthetic  in  surgical  procedures  (Kimbrough  and  others,  1985:369).  Past 
practices,  such  as  inadequate  disposal  methods,  have  led  to  a  problem  with  TCE 
contaminated  groundwater.  In  fact,  TCE  is  the  most  common  organic  groundwater 
contaminant  and  occurs  at  higher  concentrations  than  any  other  organic  pollutant 
(Marrin  and  Thompson,  1987:21).  With  its  numerous  and  widely  varied  operations, 
the  Air  Force  has  found  that  it  has  its  own  problem  with  TCE  contamination  (Gross 
and  Termaath,  1985:119).  One  of  the  consequences  of  this  contamination  is  the 
potential  for  the  contamination  of  food  crops,  whether  in  commercial  farms  or 
private  garden  plots,  and  the  subsequent  exposure  of  humans.  At  McClellan  Air 
Force  Base,  TCE  originating  from  Air  Force  operations  has  been  detected  in  off-base 
private  wells  that  are  being  used  to  irrigate  garden  plots  (ATSDR,  1994:33).  Hill 
Air  Force  Base  has  a  similar  problem  with  irrigation  of  food  crops  using 
groundwater  contaminated  by  TCE  from  the  base  (Radian,  1994:8-2). 

The  current  approach  to  identifying  crop  exposure  problems  focuses 
primarily  on  sampling  food  products  for  the  presence  of  chemicals.  Very  little  has 
been  done  to  examine  the  ways  that  organic  chemicals  are  transported  to  and  taken 
up  by  the  crops.  If  more  was  known  about  the  ways  in  which  TCE  becomes 
available  to  the  plants  for  uptake,  problems  could  be  more  readily  identified  or 


anticipated  and  sampling  and  remediation  efforts  could  be  directed  more  efficiently. 

The  goal  of  this  thesis  is  to  develop  a  methodology  for  assessing  the  influence 
of  environmental  conditions  on  the  potential  exposure  of  a  food  crop  to 
trichloroethylene  from  a  contaminated  aquifer.  Specificly,  the  methodology  will 
address  a  mature  com  crop.  This  methodology  will  take  the  form  of  a  series  of 
environmental  compartments  which  will  be  modeled  to  evaluate  the  maximum 
potential  eiqposure.  In  order  to  assess  the  usefulness  of  the  information  obtained, 
these  models  will  be  applied  to  the  example  of  Hill  AFB. 

Figure  1  shows  the  physical  cnviromnent  to  be  assessed.  The  first  step  is  to 
evaluate  the  pathways  TCE  may  take  from  the  aquifer  to  the  crop  and  the  significant 
fates  of  TCE  along  those  pathways.  The  TCE  is  assumed  to  start  in  the 
contaminated  aquifer.  It  may  be  transported  in  vapor  form  within  the  air  space  of  the 
soil  to  the  area  beneath  the  plants.  From  there,  it  potentially  continues  on  to  the  air 
of  the  crop  canopy  to  form  an  air  ejqposure.  Another  pathway  occurs  when  the 
contaminated  water  is  pumped  to  the  surface  and  used  for  irrigation.  One  of  two 
irrigation  methods  might  be  used,  spray  irrigation  or  flood  irrigation.  Once  it 
reaches  the  crop,  the  TCE  may  vaporize  into  the  air  of  the  crop  canopy  or  it  may 
leach  into  the  soil  surrounding  the  plant  roots.  The  next  step  of  the  methodology 
will  be  to  estimate  these  exposures  in  the  air  and  soil. 

After  the  development  of  the  methodology,  results  associated  with  Hill  AFB 
will  be  examined.  Where  actual  data  is  not  available,  the  conditions  at  Hill  AFB 
will  be  approximated  as  nearly  as  possible.  By  examining  the  results  from  the 
models,  the  type  of  information  to  be  gained,  as  well  as  its  usefulness,  will  be 
ascertained. 

This  methodology  will  be  a  first  step  in  better  understanding  situations  A^iiere 
crop  exposure  to  aquifer  contaminants  is  a  problem.  Based  on  the  information 
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gained,  sampling  regimes  may  be  developed  more  efficiently.  The  methodology  may 
also  reveal  those  situations  where  the  exposure  may  be  considered  negligible. 

The  series  of  models  may  also  provide  information  by  which  the  risk  to 
humans  may  be  decreased.  By  assessing  the  significant  factors  which  lead  to  crop 
exposure,  the  models  will  provide  insight  into  ways  in  which  potential  problems  can 
be  corrected.  By  contributing  information  and  understanding,  this  methodology  will 
be  a  potentially  valuable  tool  in  identifying  and  decreasing  the  overall  risk  to  human 
health. 
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II.  Literature  Review 


Trichloroethylene 

Chemical  Properties.  The  chemical  properties  of  TCE  determine  how  it 
behaves  in  the  environment.  Some  of  the  more  important  parameters  used  to 
describe  the  behavior  of  organic  chemicals  are  the  octanol/water  partition  coefficient, 
the  organic  carbon  partition  coefficient,  vapor  pressure,  solubility,  and  Henry's  Law 
constant  The  values  of  these  properties  and  others  of  inq>ortance  are  listed  in  Table 
1 .  These  {H'operties  not  only  serve  as  relative  measures  of  a  certain  behavior  but  are 
also  used  to  analytically  describe  and  predict  particular  behaviors  of  interest 

OctanolAVater  Partition  Coefficient  The  octanol/water  partition 
coefficient  (K^vv)  ^  measurement  of  the  relative  solubility  of  a  chemical  in  octanol 
versus  water.  A  large  KQ^.generally  indicates  a  tendency  to  accumulate  in  soil, 
sediments,  and  the  lipid  phase  of  biological  systems.  On  the  other  hand,  a  chemical 
with  small  partition  into  the  air  or  water  (ATSDR,  1992:6-12). 

Trichloroethylene  has  a  relatively  low  Kqvv  (log  K^vv  =2.29XBourg  and  others, 
1992:359).  Due  to  this  low  bioaccumulation  in  the  food  chain  is  not  thought 
to  be  a  significant  concern  (Callahan  and  others,  1979:52-4).  This  small  coefficient 
value  also  reflects  the  fact  that  TCE  is  weakly  hydrophobic  (Bourg  and  others, 
1992:359). 

Organic  Carbon  Partition  Coefficient.  The  organic  carbon  partition 
coefficient  (Kq^)  measures  the  tendency  of  a  chemical  to  be  adsorbed  by  soil  and 
sediment  (ATSDR,  1992:6-12).  This  is  a  chemical  specific  parameter  (ATSDR, 
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1992:6-1 1)  but  it  may  change  dramatically  with  the  particular  properties  of  the 
individual  soil  (Domenico  and  Schwartz,  1990:443).  Because  of  the  strong 
dependence  of  K^c  tbe  individual  soil  being  dealt  ?  value  is  not  given  here.  There  is 
a  strong  correlation  between  K^^andK^^  because  partition  between  organic 
carbon  and  water  is  very  similar  to  partitioning  between  octanol  and  water 
(Domenico  and  Schwartz,  1990:443). 

Vapor  Pressure.  The  vapor  pressure  of  TCE  is  important  because  it 
indicates  the  tendency  of  the  chemical  to  volatilize  from  the  liquid  phase  to  the 
gaseous  phase.  Chemicals  with  a  high  vapor  pressure  have  the  ability  to  move  from 
liquid  phase  in  the  soil  to  the  soil  air  phase  where  they  may  be  quite  mobile.  The 
v24X)r-phase  may  be  taken  up  by  the  plant  below  ground  by  the  roots  or  it  may 
volatilize  to  the  ambient  air  and  be  taken  up  by  the  above  ground  portions  of  the 
plant  (Ryan  and  others,  1988:2305). 

The  vapor  pressure  is  the  pressure  whidi  exists  in  the  gas  phase  when  the  gas 
and  liquid  phases  of  the  chemical  are  in  equilibrium  (Weast,  1985:F105).  As 
temperature  increases  the  vapor  pressure  will  increase.  On  a  molecular  level,  this  is 
because  that  increases  in  temperature  increase  the  average  kinetic  energy  of  the 
molecules  and  subsequently  increase  the  number  of  molecules  with  sufficient  kinetic 
energy  to  escape  the  liquid  into  the  gaseous  phase  (Bodner  and  Pardue,  1989:G-1 5). 
The  vapor  pressure  of  TCE  is  74  mm  Hg  at  20®C  (Dilling,  1977:407). 

The  solubility  of  a  chemical  is  a  measurement  of  how  much  of  the  substance 
will  dissolve  into  solution  with  a  certain  solvent  and  is  normally  described  with 
water  as  the  solvent.  Solubility  of  non-ionic  chemicals  is  mostly  determined  by  the 
relative  polarity  of  the  solute  and  solvent.  Like  tends  to  dissolve  like.  For  example, 
weakly  polar  solutes  are  most  effectively  dissolved  by  weakly  polar  solvents 
(Morrison  and  Boyd,  1992:32).  The  solubility  of  TCE  in  water  has  been  measured 
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as  1 .366  g/L  at  20®C  (ATSDR,  1989:33).  This  is  considered  a  low  solubility.  It 
may  be  implied  that  this  is  because  TCE  is  considered  weakly  polar  or  non  polar  and 
water  is  highly  polar. 

Henry's  Law  Coefficient.  The  Henry's  Law  coefficient  (Kjj)  is  a 
constant  which  measures  the  tendency  of  a  chemical  to  volatilize  from  a  solution. 

The  molecular  weight,  vapor  pressure,  and  solubility  all  come  into  play  in 
determining  how  strongly  the  chemical  is  driven  into  solution  relative  to  the  gaseous 
phase,  and  they  all  help  determine  the  value  of  this  constant  (ATSDR,  1992:6-9).  If 
something  happens  to  decrease  the  solubility  of  a  chemical,  such  as  a  change  in 
temperature,  the  change  is  reflected  in  this  coefficient.  In  a  similar  manner,  if  the 
v^)or  pressure  of  the  chemical  changes,  the  tendency  to  volatilize  from  the  solution 
changes.  Once  again,  the  Henry's  Law  coefficient  reflects  the  new  relationships 
(Bodner  and  Pardue,  1989:597). 

In  contrast  to  vapor  pressure,  which  measures  the  pressure  of  the  liquid,  this 
coefficient  indicates  the  partitioning  concentrations  between  the  liquid  and  gas  phase 
at  equilibrium.  Analytically,  Henry's  Law  can  be  expressed  as: 

Cg=KHCL 

where  Cg  and  are  the  concentrations  in  the  gas  and  liquid  phases  respectively.  A 

substance  with  a  higher  Henry's  Law  constant  moves  to  the  gaseous  phase.  The 
Henry's  Law  constant  for  TCE  is  1.18  (kPa  m^/mol  at  25®C)  For  volatilization  to 
be  significant,  an  organic  solvent  should  have  a  vapor  pressure  of  1 30  Pa  or  greater 
and  a  Henry's  Law  constant  of  50  Pa  m^/mol,  as  TCE  does  (Bourg  and  others, 
1992:360). 

Transformation.  The  loss  of  trichloroethylene  to  degradation,  reaction,  or 
other  transformation  is  important  to  knowing  a  compound's  persistence  in  the 
environment.  Whether  or  not  it  degrades  under  certain  conditions  determines 


7 


whether  or  not  it  makes  it  to  the  receptor,  in  this  case  the  com  plant,  and,  if  it 
remains  in  the  vicinity  of  the  receptor,  whether  or  not  it  accumulates.  In  the  problem 
under  consideration. 


Parameter 

Standard  Svmbol 

Value 

Source 

Octanol  Water 
Partition  Coefficient 

Kow 

log  Kqyv  “  2.29 

(Bourg  and  others, 
1992:359) 

Vapor  Pressure 

9.8  kPa 
at20'»C 

(Dilling, 

1977:407) 

Aqueous  Solubility 

1.366  g/L  at  20‘’C 

(ATSDR, 

1989:33) 

Henry's  Law  Constant 

Kh 

1.18kPam3/mol 

at25®C 

(Bourg  and  others, 
1992:360) 

Adsorption 

Equilibrium  Constant 
for  TCE 

Kd 

11.7mL/g 

(Pavlosthis  and 
offiers,  1991:277) 

Gaseous  Molecular 
Difiusion  Coefficient 

D 

7030  cm2/day 

(Peterson  and 
others,  1988:576) 

Molecular  Weight 
TCE 

M 

130  g/mol 

Table  1:  Summary  ofchemical  properties  of  TCE. 


there  are  several  environments  in  ^^ch  TCE  is  being  examined  and  may  be 
degraded  or  react  to  form  other  chemicals.  The  air,  unsaturated  soil,  and  the 
aqueous  environments  are  all  separate  compartments  in  this  model  where  this  may 
luqjpen. 
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Reaction  with  hydroxyl  radicals  is  the  primary  loss  pathway  for 
trichloroethylene  in  air  (USEPA,  1985;  1-1).  Hydroxyl  radicals  form  in  the 
atmosphere  in  the  following  reaction  involving  ozone  and  water  vapor: 

O3  +  hv  0(D')  +  O2 
0(D')  +  H2O  ->  20H 

where  hv  represents  ultraviolet  light  and  CXD')  is  excited  oxygen.  By  the  above 
equations,  it  is  apparent  that  the  hydroxyl  radicals  are  formed  in  the  presence  of 
sunlight  (Godish,  1991 .50)  and  the  atmospheric  lifetime  of  the  hydroxyl  radicals  is 
on  the  order  of  microseconds  (Hauser,  1994). 

This  is  important  m  that  the  crop  canopy,  the  air  environment  of  concern,  is 
assumed  to  not  allow  a  substantial  amount  of  sunli^t  in.  It  is  doubtful,  considering 
the  atmospheric  lifetime  of  the  hydroxyl  radicals,  that  any  significant  quantity  of 
hydroxyl  radicals  would  persist  long  enough  to  be  transported  to  portions  of  the 
canopy  not  radiated  by  the  sun.  This  is  perh^  a  conservative  assunq>tion  because 
some  sunlight  would  reach  some  levels  of  the  crop  below  the  initial  canopy  and  some 
reactions  would  therefore  take  place.  In  the  troposphere,  half  life  of  TCE  has  been 
estimated  from  54  hours  to  2  weeks(USEPA,  1985:3-6).  No  studies  have  been 
found  which  are  iq)plicable  to  an  air  environment  which  lacks  con:q>lete  ultraviolet 
light. 

The  second  environment,  the  unsaturated  soil  environment,  is  even  more 
difficult  to  quantify.  A  review  of  the  litCTature  has  not  uncovered  direa  and 
quantifiable  indications  degradation  or  other  loss  of  TCE  in  the  soil  environment 
This  may  be  due  to  the  inability  to  quantify  an  actual  source  of  the  TCE  with 
sufficient  accuracy  as  well  as  difficulty  in  quantifying  loss  pathways,  such  as 
partitioning  into  the  soil.  In  essence,  the  necessary  mass  balances  have  not  been 
established.  Various  studies  that  have  attempted  to  estimate  degradation  in  soil 
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systems  have  postulated  biodegradation  as  the  most  likely  transformation  process  but 
have  estimated  it  to  be  slow  if  present  at  all  (ATSDR,  1989:9 1 ). 

Flood  irrigation  often  uses  open  trenches  to  transport  the  water  to  the  plant. 

At  that  time  the  TCE  is  in  a  purely  aqueous  environment.  There  are  no  loss 
pathways  which,  in  an  aerated  aqueous  environment,  are  significant  compared  to 
volatilization  (ATSDR,  1989:91). 

Soil  Gas  Pathway 

Sorption.  Sorption  is  a  term  ^^ch  refers  in  general  to  a  situation  where  the 
soil  matter  in  some  way  binds  up  the  molecules  of  a  contaminant.  The  two 
mechanisms  which  normally  categorize  sorption  are  adsorption  and  absorption.  In 
adsorption,  the  contaminant  physically  adheres  to  the  exterior  surface  of  the  sorbing 
material.  In  absorption,  the  contaminant  is  taken  into  the  interior  of  the  soil  matter 
in  a  manner  similar  to  the  way  a  sponge  would  soak  up  water.  In  both  cases  the 
contaminant  is  no  longer  mobile,  and  its  activity  in  the  environment  is  inhibited. 
Sorbed  chemical  is  also  not  available  for  uptake  by  the  plant  (Boesten,  1993:405). 
Hence,  the  sorption  of  TCE  to  soil  matter  has  a  primary  role  in  determining  the 
activity  of  the  chemical  in  the  soil  environment  (Chiou,  1990: 1 52). 

Soil  has  two  solid  phases  which  play  a  role  in  sorption,  the  mineral  portion 
and  the  organic  portion.  The  mineral  cont^t  of  the  soil  acts  as  a  conventional  solid 
adsorbent,  and  the  soil  organic  matter  acts  as  a  partition  medium,  or  absorbent 
(Chiou,  1990: 1 1 8).  Only  a  finite  amount  of  water  and  contaminant  can  adsorb  to 
the  surface  of  a  given  amount  of  mineral.  Water  is  more  polar  than  TCE,  and  this 
increased  polarity  gives  it  a  stronger  tendency  to  be  adsorbed  (Chiou  and  others, 

1 985 : 1 1 98).  In  the  resulting  competition,  the  moisture  in  the  soil  decreases  the 
amount  of  TCE  which  can  be  adsorbed. 
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The  moisture  content  of  the  soil  is  therefore  an  important  factor  in  the  mineral 
sorption  activity  of  TCE,  and  mineral  sorption  becomes  a  less  dominant  mechanism 
as  moisture  content  increases  (Peterson  and  others,  1988:571).  Because  of  this 
competition  with  the  water,  sorption  of  TCE  in  an  unsaturated  environment  is  quite 
different  from  sorption  of  TCE  in  a  saturated  environment  (Peterson  and  others, 
1988:571).  Water  competes  so  strongly  for  adsorption  to  the  minerals  that  it  can 
actually  displace  TCE  from  the  mineral,  and  an  increase  in  moisture  content  can 
cause  a  release  of  TCE  (Spencer  and  others,  1982:20).  The  water  content  necessary 
to  sustain  plants,  known  as  the  wilting  point  is  sufficient  to  displace  the  water  from 
the  mineral  content  and  make  partitioning  to  the  soil  organic  content  the  dominant 
sorption  mechanism  (Chiou,  1990:148).  The  wilting  point  is  the  moisture  content 
below  which  plants  wilt  and  do  not  recover  (Hartman  and  others,  1988:666).  Thus, 
mineral  sorption  is  not  a  controlling  mechanism  in  a  healthy  agricultural 
environment 

TCE  partitions  into  the  soil  organic  content  by  forces  typical  to  a  substance 
going  into  solution,  such  as  Van  der  Waals  forces  (Chiou,  1990: 1 52).  Van  der 
Waals  forces  are  the  electrostatic  interactions  between  molecules.  Sorption  by  soil 
organics  increases  linearly  with  the  TCE  concentration.  This  is  in  contrast  to 
sorption  by  the  mineral  portion,  which  is  very  nonlinear  (Chiou  and  others, 

1985 : 1200).  As  a  result,  sorption  ranges  from  being  highly  linear  to  highly 
nonlinear  depending  on  the  relative  dominance  of  mineral  and  organic  sorption.  The 
total  adsorption  activity  of  the  soil/TCE  system  is  affected  by  such  things  as  the 
properties  of  TCE,  TCE  concentration,  soil  moisture  content,  and  the  soil  properties, 
such  as  amount  of  mineral  or  organic  matter  (Spencer  and  others,  1982: 17). 

Flux  Rates.  Flux  rate  refers  to  the  vapor  flux  from  the  contaminated  aquifer 
to  the  crop  environment  through  the  soil  atmosphere,  or  soil  air  space.  This  flux 
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continues  to  the  canopy  and  on  to  the  atmospheric  sink.  There  are  several  important 
pieces  in  this  flow  process.  The  first  is  the  volatilization  of  the  TCE  from  the 
contaminated  aquifer  to  the  soil  air.  The  second  is  the  transport  of  the  vapor  through 
the  porous  soil  to  the  soil/air  interface.  From  the  soil/air  interface  the  TCE  difruses 
into  the  crop  canopy  and  thus  becomes  the  first  of  the  two  exposure  pathways  being 
examined. 

The  initial  step  in  this  pathway  is  the  volatilization  of  the  TCE  from  the 
contaminated  aquifer  to  the  soil  atmosphere.  One  determinant  in  this  process  is  the 
amount  of  contaminant  in  the  capillary  fringe  (Marrin,  1988:741).  The  capillary 
fringe  is  a  saturated  zone  over  the  aquifer.  Water  is  drawn  up  into  this  zone  by  the 
capillary  attraction  of  the  small  pore  spaces  in  the  soil  (Fetter,  1988:91). 

Capillary  attraction  is  the  slight  negative  pressure  which  develops  in  the  pore 
spaces  in  contact  with  the  water  (Domenico  and  Schwartz,  1990:89)  due  to  the 
molecular  attraction  between  the  soil  particles  and  the  water  (Fetter,  1 988:9 1 ).  At 
the  capillary  fringe  the  chemical  can  come  into  contact  with  the  air  and  volatilize.  It 
has  been  hypothesized  that  the  chemical  may  move  up  through  the  capillary  fringe 
by  a  combination  of  dispersion  and  molecular  difiusion.  In  addition,  water  table 
fluctuations  may  greatly  increase  the  amount  of  volatilization  (Marrin,  1988:741). 
TCE  is  more  dense  than  water  and,  tending  to  sink,  may  be  found  in  higher 
concentrations  progressively  lower  in  the  aquifer  (USEPA,  1992: 1).  A  decrease  in 
the  water  level  may  bring  the  capillary  fringe  into  contact  with  this  higher 
concentration  and  allow  greatly  increased  volatilization. 

After  vaporizing  from  the  aquifer,  the  vapor  phase  TCE  is  in  the  soil  gas 
phase  in  the  unsaturated  soil  environment  In  general,  movement  of  vapor  in  the 
unsaturated  zone  is  attributed  to  molecular  difiusion  (Bishop  and  others,  1990:257). 
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Molecular  diffusion  is  the  spreading  due  to  random  molecular  motion  and  collisions 
(Jury,  1986: 135).  Flux  of  a  gas  due  to  diffusion  can  be  expressed  as: 

F=  -D  dC/dz 

and  is  this  is  known  as  Fick's  First  Law.  F  is  the  mass  flux  per  unit  area  per  unit 
time.  D  is  the  diffusion  coefficient  in  units  of  area  per  unit  time.  dC/dz  is  the 
concentration  gradient  of  the  vapor.  The  negative  sign  indicates  that  the  flow  is  from 
greater  to  lesser  concentrations  (Fetter,  1988:390).  Transport  by  diffusion  implies 
that  the  concentration  in  the  soil  gas  would  increase  with  depth  and  with  proximity  to 
the  larger  concentrations  of  the  source.  This  has  been  observed  in  at  least  one  case 
study  (Marrin  and  Kerfoot,  1988:741). 

Whether  by  directly  influencing  the  diffusion  coefficient  of  the  chemical  or  by 
in^acting  the  transport  in  some  other  way,  many  elements  can  affect  the  rate  of 
transport  in  the  soil  gas  phase,  such  as: 

-the  porosity  of  the  soil 

-the  Henry's  Law  properties  of  TCE 

-the  adsorption  to  the  soil 

-the  water  content  of  the  soil 

-the  evaporation  of  water  content 

-the  variations  in  the  soil  composition 

All  of  these  factors  combine  to  determine  the  amount  of  chemical  available  and  the 
rate  at  which  it  is  transported  through  the  soil. 

The  porosity  of  the  soil  is  defined  as  the  amount  of  air  space  in  the  soil  matrix 
and  is  a  property  which  is  characteristic  of  the  soil  type.  It  is  deflned  mathematically 
as 

(j)=  100  VyA^ 
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where  (|>  is  the  porosity  in  percentage,  Vy  is  the  volume  of  the  void  space  in  a  unit 
volume  of  soil,  and  V  is  the  unit  volume  of  soil  with  both  the  solids  and  the  void 
space  (Fetter,  1988:63-64).  This  same  property  is  often  referred  to  as  bulk  density, 
which  is  just  the  antithesis.  As  bulk  density  increases,  porosity  decreases  and  vice 
versa.  Porosity  affects  the  rate  at  which  vapor  chemical  is  transported  by 
determining  the  cross  sectional  area  available  for  diffusion.  It  also  increases  or 
decreases  the  path  length  that  the  diffusing  molecules  must  travel  by  determining  the 
amount  of  solids  per  unit  volume  that  the  vapor  must  go  around.  Porosity  can  be 

incorporated  into  Pick's  Law  by  a  tortuosity  factor.  This  is  expressed  as; 

Jy  =  -Ti  D  5Cg  /  5z 

where  Jy  is  the  vapor  ffux.  t]  is  a  tortuosity  factor  which  accounts  for  changes  in 
cross  sectional  area  and  path  length  due  to  changes  in  porosity. 

Chemicals  in  the  soil  environment  are  governed  by  Henry's  Law  and  its 
implications  for  volatilization  rates  (Jury(b),  1986: 1 59).  The  Henry's  Law 
properties  of  the  chemical  vary  with  enviroimiental  factors,  such  as  soil  temperature. 
This  may  be  most  important  at  the  capillary  fringe  \^ere  the  chemical  is  initially 
vaporizing. 

Adsorption  is  important  in  that  it  decreases  the  amount  of  chemical  which  can 
be  transported  to  the  surface  and  subsequently  to  the  canopy  environment  (Jury(b), 
1986: 167).  It  forms  an  additional  manner  in  which  TCE  may  be  lost  and  not  expose 
the  crop. 

Water  content  has  several  potential  effects  and  these  effects  may  be 
contradictory.  Soil  moisture  content  will  tend  to  decrease  the  transport  of  TCE  to  the 
surface.  Water  content  fills  the  pore  spaces  and  decreases  the  amount  of  air  space 
within  the  soil  available  for  diffusion  (Jury,  1 986;  1 36).  If  water  is  evaporating 
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within  the  soil  matrix  it  may  carry  TCE  with  it  by  mass  flow  and  cause  an  increase 
in  the  transport  of  the  contaminant.  This  effect  may  be  less  for  TCE  than  for  other 
chemicals  due  to  its  relative  insolubility  (Spencer  and  others,  1982:21).  Water 
evaporation  at  the  surface  and  increased  use  of  water  by  plants  may  cause  an 
increase  in  the  gradient  of  water  vapor  and  a  subsequent  increase  in  ditfusion  of 
water  vapor  (Taylor  and  others,  1976:629).  The  reasoning  is  apparent  by 
observation  of  Pick's  Law  as  diflusion  increases  proportionally  with  the 
concentration  gradient.  This  increase  in  the  gradient  due  to  evaporation  is  known  as 
the  "wick  effect"  (Spencer  and  others,  1982:22). 

Correlation  with  Groundwater  Concentration.  While  a  significant  amount  of 
research  has  been  accumulated  on  flow  in  saturated  soils,  relatively  little  is  known 
about  contaminant  flow  in  soil  gas  and  other  unsaturated  systems  (Peterson  and 
others,  1988:571).  The  variations  in  subsurface  conditions,  many  of  which  are 
difficult  to  characterize,  make  quantitative  relationships  between  a  groundwater 
contamination  and  a  soil  gas  contamination  difficult  to  establish  (Bishop  and  others, 
1990:257).  A  review  of  the  current  literature  has  uncovered  no  predictive  models  of 
this  environment.  Several  efforts  have  been  made  to  statistically  correlate 
groundwater  contamination  to  soil  gas  contamination,  and  only  order  of  magnitude 
relationships  have  been  established  (Bishop  and  others,  1990:257). 

Soil-Air  Interface.  Once  the  TCE  reaches  the  surface,  volatilization  is 
dependent  on  the  vapor  pressure  of  the  chemical  and  its  rate  of  movement  away  from 
the  surface  (Spencer  and  others,  1988:504)  At  the  soil  surface,  the  atmosphere 
generally  has  some  turbulence  wliich  remo  ves  the  TCE  rapidly  but,  adjacent  to  the 
surface,  there  is  a  stagnant  boundary  through  which  the  vapor  must  move  by 
molecular  diffusion  (Jury,  1 986b;  161).  Rate  of  diffusion  through  this  stagnant  layer 
is  an  important  and  possibly  limiting  step  in  the  volatilization  to  the  crop  canopy. 
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Precipitation  tends  to  decrease  volatilization  because  the  water  tends  to  move 
chemicals  away  from  the  surface  by  mass  flow  as  it  penetrates  down  through  the  soil. 
Irrigation  with  uncontaminated  water  has  an  effect  similar  to  that  of  precipitation 
(Jury,  1986b:  169).  Irrigation  and  precipitation  are  important  limiting  factors  on  the 
volatilization  of  soil  gas. 

Numerous  emperical  models  exist  for  predicting  the  volatilization  of  soil 
incorporated  chemicals.  Much  of  this  work  was  accomplished  for  the  purpose  of 
characterizing  the  behavior  of  pesticides  in  the  soil  environment.  Unfortunately, 
many  of  the  parameters  used  in  these  model  have  to  be  determined  experimentally 
and  are  not  available  for  TCE.  The  models  were  developed  and  tested  for 
agricultural  chemicals,  and  TCE,  lacking  agricultural  use,  was  not  included  in  the 
studies. 

Volatilization  Related  Fates 

Volatilization  From  Ponded  Water  During  irrigation  with  contaminated 
groundwater,  a  certain  amount  of  TCE  will  be  residual  in  the  water  after  it  reaches 
the  crop  canopy,  and  the  water  would  be  expected  to  pond  for  a  period  of  time. 
During  that  time,  the  TCE  will  be  transf^ed  into  the  air  of  the  crop  canopy. 
Volatilization  of  TCE  from  a  body  of  water  can  be  described  by  a  two  layer  model  as 
shown  in  Figure  2. 

The  liquid  film  and  gas  film  parts  form  the  interface.  The  films  are  assumed 
to  be  stagnant,  and  therefore  molecular  difliision  is  the  primary  mechanism  for 
movement  of  the  chemicals  to  or  away  from  the  interface  (Liss  and  others, 

1974: 181).  These  diffusion  coefficients  have  been  evaluated  for  TCE.  The  bulk  gas 
phase  and  bulk  liquid  phase  are  part  of  the  main  body  of  atmosphere  or  water. 
Movemert  of  the  chemical  in  this  phase  is  controlled  by  turbulence  and  molecular 
difhision.  Where  turbulence  is  significant,  it  would  be  the  dominant  mechanism.  In 
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Turbulent  and  Molecular  Transfer 


Bulk  air  phase 


Gas  film 

^  Molecular  Transfer 

Liquid  film 

^  Molecular  Trar«sfer 

Interface 


Turbulent  and  Molecular  Transfer  t  Bulk  liquid 

phase 

Figure  2:  Two-layer  model  of  a  gas-llquid  Interface  (LIss  and  others,  1974:181) 


the  case  of  water  which  has  ponded  under  the  crop  canopy,  turbulence  will  probably 
be  negligible.  The  volatilization  of  the  chemical  requires  that  it  first  be  transported  to 
the  liquid  film  layer,  and  it  must  then  difiuse  through  the  interface.  At  the  interface 
the  volatilization  is  controlled  by  Henry's  Law  (Spencer  and  others,  1988:505).  The 
two  rate  controlling  mechanisms  are  difiiision  in  the  films  and  volatilization  at  the 
interface.  The  slower  of  the  two  mechanisms  will  control  the  rate  at  which  the 
contaminate  enters  the  atmosphere  (Mackay  and  Wolkoff,  1973 :6 1 1 ).  This  does 
assume  that  the  transfer  of  the  chemical  is  not  rapid  enou^  to  cause  depletion  of  the 
chemical  near  the  liquid  film.  Contaminates  that  have  a  high  vsqior  pressure  tend  to 
be  limited  in  their  volatilization  by  the  difiiision  through  the  liquid  film  as  the 
transfer  at  the  interface  occurs  more  readily  (Mackay  and  Leinonen,  1 975 : 1 1 78). 
One  of  the  unexpected  aspects  of  the  volatilization  of  a  dilute  solute  such  as  TCE  in 
water  is  that  the  chemical  will  evaporate  more  rapidly  than  the  water  solvent  even 
though  the  concentrations  are  within  the  solubility  limits  of  the  chemical  (Dilling 
and  others,  1975:833). 

Another,  slightly  more  complex  situation  where  volatilization  occurs  in  the 
problem  under  review  would  be  the  flow  in  the  open  ditch  associated  with  transport 
to  the  crop  during  flood  irrigation.  Volatilization  has  been  studied  empirically 
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under  many  specific  aqueous  environments.  They  are  distinguished  by  the  factors 
which  control  the  rate  at  which  the  TCE,  or  other  solute,  vaporizes.  Such  factors 
may  include  the  following  (Chiou,  1990;  Liss  and  others,  1974;  Barrio-Lage  and 
others,  1987;  Mackay  and  Leinonen,  1975): 


Depth:  The  depth  of  the  solvent  determines  the  maximum  distance 
which  the  solute  will  have  to  be  transported  to  reach  the  liquid 
film. 

Surface  Area:  The  surface  area  of  the  solvent  determines  the  size 
of  the  interface  available  to  take  part  in  the  transfer. 

Air  Flow:  If  the  air  flow  in  the  bulk  air  phase  is  not  sufficient  to 
prevent  a  significant  vapor  pressure  from  developing,  the 
vaporization  will  be  limited  according  to  Henry's  Law. 

Turbulence:  Turbulence  in  the  bulk  liquid  prevents  depletion  of 
the  chemical  near  the  liquid  film  from  limiting  the  rate  of 
volatilizatioa 

Difiusion  Coefficient:  The  rate  at  which  a  chemical  difiuses 
through  the  air  and  gas  film  layers  can  control  the  rate  at  which  the 
chemical  is  transferred  to  the  air. 

Chemical  Properties:  The  specific  properties  of  the  chemical,  such 
as  its  Henry's  Law  constant  and  vapor  pressure  are,  as  stated 
previously,  inqxntant  in  determining  the  rate  of  vaporization. 


Because  of  these  distinguishing  elements,  studies  may,  for  example,  measure  one 
rate  of  volatilization  for  a  solute  in  a  vessel  stirred  at  a  certain  rate  and  a  very 
different  rate  for  volatilization  from  a  lake. 

Aeration  and  Volatilization.  Spray  irrigation  has  the  unique  aspect  of 
propeling  the  water  through  the  air  in  droplet  form.  This  act  provides  a  unique 
opportunity  for  the  TCE  to  vaporize  and  to  be  lost  to  the  troposphere.  Nearly  all  of 
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the  information  available  about  the  behavior  of  TCE  under  these  circumstances  was 
collected  for  the  purpose  of  studying  exposure  in  an  indoor  environment  from 
residential  use  of  contaminated  water,  especially  in  a  shower. 

One  of  the  most  important  considerations  in  the  volatilization  from  drops  is 
the  drop  size.  Drop  size  is  important  because  it  controls  the  surface  area  of  the 
interface  between  the  liquid  and  the  air  phases  (Giardino  and  others,  1992: 1602). 
The  amount  of  volatilization  is  also  impacted  by  the  chemical  concentration  in  the 
gas  form  surrounding  the  drop,  as  might  be  expected  from  Henry's  law. 

The  residence  time  of  the  drop  in  the  air  controls  the  amount  of  time  the 
chemical  has  to  vaporize  before  the  drop  reaches  the  end  of  its  trajectory  and  changes 
phase  (Giardino  and  others,  1992: 1603).  The  most  complex  element  may  be  the 
mass  transfer  coefficient.  This  is  analogous  to  the  transfer  through  the  liquid  film 
phase  of  the  simple  air/water  interface  in  Figure  2  above.  It  is  often  treated  as  a 
resistance  to  volatilization. 

There  are  several  theories  which  attempt  to  explain  the  conditions  of  the 
drop  and  how  the  chemical  interfaces  with  the  surrounding  air.  These  theories  range 
from  turbulent  conditions  within  the  drop  to  a  drop  ufrich  oscillates  and  changes 
shape  during  oscillation.  Drop  size  is  often  a  parameter  in  theoretical  attempts  to 
estimate  the  mass  transfer  coefficient  (Altwicker  and  others,  1988:329).  It  is 
noteworthy  that  for  chemicals  with  a  high  Henry's  Law  constant,  such  as  TCE, 
McKone  and  Knezovich  (1991  ;836)  found  the  transfer  efficiency  to  be  a  strong 
function  of  the  rate  of  difilision  in  water  and  not  a  function  of  temperature.  In 
addition,  Altwicker  (1988:33 1)  concludes  that  transfer  to  the  gas  phase  is  not 
velocity  dependent  despite  the  velocity  dependence  of  most  of  the  theories  which 
attempt  to  explain  the  transfer. 
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Pumping  Losses 

Prior  to  being  used  for  irrigation  purposes,  water  must  be  raised  from  the 
aquifer.  This  is  typically  done  with  one  of  a  wide  variety  of  pumps.  A  seemingly 
limitless  number  of  parameters  have  been  studied  as  to  how  they  inqiact  the 
concentration  of  TCE  which  is  delivered  by  the  pumping  system.  These  include 
everything  from  pumping  rate  to  the  type  of  piping  used.  The  original  intent  of  these 
studies  was  to  study  the  representativeness  of  samples  taken  for  monitoring  purposes, 
and  the  equipment  being  evaluated  is  quite  different  than  that  which  would  be  used 
for  residential  or  commercial  water  needs. 

Many  pumps  lift  the  water  by  creating  a  negative  suction.  From  Henry's 
Law,  this  might  be  expected  to  cause  a  volatile  compound  to  enter  the  air  phase. 
Pearsall  (1987:64)  found  that  this  was  not  the  case  for  lift  heights  of  up  to  24.5  feet. 
Ho  (1983:586)  concluded  that  this  pressure  drop  did  cause  a  loss  of  TCE  at  lift 
heights  greater  than  16  feet  and  that  it  increased  with  increasing  lift  height 

Pumping  rate  is  often  cited  as  an  important  factor  in  determining  the  amount 
of  TCE  in  the  water  which  is  delivered.  Ho  (1983:586)  found  that  increased 
pumping  rates  lead  to  significant  losses  of  TCE  but  did  not  describe  a  potential 
causality. 

The  materials  used  in  the  pumping  system  may  also  be  important  in 
determining  the  amount  of  TCE  which  is  delivered.  For  example,  metals  strongly 
adsorb  organic  compounds  and  act  as  a  catalyst  in  some  reactions  (Pettyjohn  and 
others,  198 1 : 1 8 1).  Many  plastics  will  adsorb  organics  as  well  (Pettyjohn  and 
others,  1981:181).  However,  these  losses  have  not  been  quantified  well  enough  to 
include  them  in  this  model. 

A  review  of  the  available  literature  has  not  produced  quantifiable  losses 
which  may  be  extrapolated  to  a  typical  irrigation  scenario.  This  may  be  an 
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important  blind  spot  in  the  current  state  of  the  research.  The  result  of  the  handling  of 
TCE  by  pumping  and  irrigation  equipment,  such  as  evaporation  or  catalyzed 
reactions,  could  have  significant  implications  for  assessments  of  risk  to  human 
health. 

Behavior  Under  the  Crop  Canopy 

The  behavior  of  the  contaminated  irrigation  water  once  it  is  within  the  crop 
canopy  will  take  two  important  forms.  Part  of  the  water,  along  with  the  TCE,  will 
evaporate  into  the  air  where  the  above  ground  portions  of  the  plant  will  be  exposed. 
In  a  similar  manner,  part  will  leach  into  the  soil  where  it  will  expose  the  plant  roots. 

In  the  air  of  the  crop  canopy,  the  most  applicable  study  is  that  of  Parmele 
(1972).  This  study  estimated  pesticide  fluxes  by  assuming  that  the  diffusivity,  or 
rate  of  difiusion,  of  the  chemical  is  the  same  as  the  difiusivity  of  water  vapor.  In  the 
case  of  an  adult  com  crop,  the  diffusivity  was  determined  by  distributing  the  wind 
drag  into  the  crop  canopy.  The  effect  is  an  estimate  of  the  turbulence  under  the 
vegetation  due  to  wind  above  the  vegetation.  After  vaporization,  this  diffusion  of  the 
TCE  up  through  the  crop  and  to  the  troposphere  will  be  critical  to  determining  the 
amount  of  accumulation. 

The  TCE  tainted  water  will  also  leach  into  the  soil.  In  the  literature,  infiltrate 
normally  refers  to  the  rate  at  which  ponded  water  enters  the  soil.  I.eaching  normally 
refers  to  the  downward  transport  of  water  wiiich  is  already  in  the  soil.  The  fate  of 
TCE  is  impacted  by  several  activities: 

-Uptake  by  the  crop. 

-Infiltration  into  the  soil. 

-Mass  flow  with  the  water. 

-Partitioning  into  the  soil. 
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Absorption  by  the  crop  is  a  potential  loss  pathway  for  organic  chemicals  (Ryan  and 
others,  1988:2307).  Although  the  uptake  by  the  plant  may  be  an  important  fate  of 
the  TCE,  such  issues  are  beyond  the  scope  of  this  thesis.  Assuming  no  losses  to  plant 
uptake  constitutes  a  conservative  assumption  for  a  model  which  is  estimating  the 
amount  of  TCE  in  the  soil  or  air. 

The  rate  of  infiltration  into  the  soil  determines  how  long  the  water  will  be 
ponded  on  the  surface  where  it  volatilizes  much  more  quickly.  It  is  affected  by  the 
rate  at  which  water  is  applied,  the  moisture  content  of  the  soil  at  the  start  of 
irrigation,  and  the  properties  of  the  soil.  The  infiltration  rate  is  higher  for  soils  with 
lower  initial  moisture  content  As  infiltration  continues  and  the  soil  becomes 
saturated,  the  rate  of  infiltration  becomes  slower  and  approaches  the  rate  associated 
with  saturated  soil  and  the  saturated  hydraulic  conductivity  (Letey  and  Oddson, 
1972:402).  The  hydraulic  conductivity  is  a  coefficient  of  proportionality  which 
describes  the  rate  at  which  a  fluid  will  move  through  soil  and  is  characteristic  of  the 
properties  of  the  soil  (Fetter,  1988:78).  This  discussion  assumes  that  the  rate  at 
which  the  water  is  applied  is  greater  than  the  rate  at  which  it  can  infiltrate  the  soil.  If 
that  is  not  the  case,  the  rate  at  which  the  water  is  applied  limits  the  rate  at  which  the 
water  can  infiltrate  the  soil. 

The  TCE  will  be  transported  by  mass  flow,  or  carried  with  the  water,  and  the 
movement  of  the  water  in  the  subsurface  is  important  to  the  fate  of  the  TCE  (Letey 
and  Oddson,  1972:402).  After  the  water  is  in  the  subsurface  it  will  penetrate  to  a 
certain  depth.  This  depth  is  affected  by  the  initial  moisture  content  of  the  soil,  the 
rate  of  application,  and  the  properties  of  the  soil  (Letey  and  Oddson,  1972:402). 
Loss  of  water  due  to  uptake  by  plants  and  evaporation  to  the  surface  can  create  a 
drying  zone  and  lead  to  the  wick  effect  discussed  previously  (Taylor  and  others. 
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1976:629).  This  upward  movement  of  water  can  carry  chemicals  back  toward  the 
surface  (Letey  and  Oddson.  1972:407). 

The  absorption  by  the  organic  content  of  the  soil  also  helps  to  determine  the 
fate  of  TCE.  If  the  soil  is  irrigated  with  contaminated  water,  two  factors  determine 
the  influence  of  absorption.  The  partition  coefficient  Kj  is  often  used  to  describe  the 
total  mass  of  absorbed  chemical  a  ftmction  of  the  solution  concentration,  C.  This 
can  be  written  as: 

S  =  K<iC 

where  S  is  the  mass  of  absorbed  chemical.  This  equation  assumes  that  an 
equilibrium  exists  between  the  solution  and  absorbed  concentrations.  The  rate  at 
which  the  absorption,  a,  approaches  equilibrium  is  also  important  (Oddson  and 
others,  1970:413).  The  amount  contaminant  in  solution  and  absorbed  tends  to 
remain  highest  near  the  surface  and  decrease  with  depth.  These  parameters 
determine  the  sh£q)e  this  concentration  profile.  If  uncontaminated  water  is  added 
later,  the  contaminant  will  move  through  the  soil  as  a  wave.  The  depth  of  the 
maximum  concentration  in  this  wave  is  influenced  by  IQi  and  tends  to  decrease  as 
increases.  The  a  influences  the  magnitude  of  the  maximum  concentration  in  this 


III.  Methodology' 


Exposure  Due  to  Soil  Gas 

Likely  Soil  Gas  Contamination.  The  best  method  currently  available  to 
estimate  the  soil  gas  flux  which  may  be  expected  from  a  contaminated  aquifer  is 
from  an  examination  of  the  existing  case  studies  on  TCE  contaminated  aquifers. 

The  basis  for  comparison  will  be  the  factors  which  primarily  determine  the  soil  flux, 
such  as  amount  of  contamination,  proximity  of  the  aquifer  to  the  surface,  and  the 
geology  of  the  vadose  zone. 

The  soil  conditions  at  Hill  AFB,  the  subject  of  the  case  study,  are  very 
heterogeneous  and  dijSicult  to  put  into  a  specific  category.  The  soils  ofif-base  where 
irrigation  activities  occur  are  silty  sand  with  lower  layers  of  silty  clay  (Radian, 
1994:3-25).  The  contaminated  aquifer  ranges  in  depth  from  the  soil  surface  to 
approximately  60  feet.  Samples  from  the  residential  wells  around  Hill  AFB 
measured  concentrations  up  to  1 1.5  ppb  TCE  (Radian,  1994:3-35). 

Several  studies  have  examined  the  correlation  between  a  groundwater 
contamination  and  soil  gas  concentration  for  TCE.  Lappala  and  Thompson  (1983) 
examined  a  site  in  the  San  Francisco,  California  area.  However,  the  soils  at  the  site 
range  from  fine  sands  to  coarse  gravels  and  are  not  comparable  to  those  at  the  Hill 
AFB  site  (Lappala  and  Thompson,  1983:671).  Runyon  and  Thompson  (1987) 
examined  a  TCE  contaminated  site  in  Kansas,  but  insufficient  information  on  the 
soil  makeup  is  provided  to  draw  analogies  to  the  Hill  AFB  site.  In  addition,  soil  gas 
profiles  necessary  for  obtaining  information  about  the  flux  rates  were  not  obtained. 
In  Marrin  and  Thompson's  (1984)  study,  two  sites  were  examined,  one  in  Tucson, 
Arizona  and  another  site  northern  California.  The  vadose  zone  at  the  Tucson  site 
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was  a  "dry,  gravelly-sand"  and  the  aquifer  was  approximately  1 1 5  feet  deep.  Both 
of  these  features  made  it  quite  different  from  the  situation  in  tlie  case  study.  The  site 
in  northern  California  was  a  shallow  aquifer,  approximately  30  feet  below  the 
surface,  and  the  vadose  zone  was  characterized  by  "a  soft,  wet,  sandy-silt  soil" 
(Marrin  and  Thompson,  1984: 173).  This  site  was,  among  those  studied,  the  most 
analogous  to  the  one  at  Hill  AFB.  The  flux  of  TCE  calculated  in  the  study  for  the 
northern  California  site  is  12  fig/day  per  m^  at  a  point  where  the  concentration  in  the 
groundwater  is  16  ppb  (Marrin  and  Thompson,  1984:173).  The  soil  gas  profile 
resulting  from  this  flux  is  given  in  Table  2. 


Deothfft) 

TCE  Concentration  fu 
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Air  Above 
Ground 

.001 
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.003 

3.5 

.01 

11 

.6 

14 

2 

20 

3 

26 

16 

Table  2:  Soil  gas  profile  at  the  northern  California  site  (Marrin  and 
Thompson,  1984:180) 
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Volatilization  From  Soil  Gas.  From  the  information  in  the  study  it  is 
possible  to  assume  that  the  transfer  from  the  soil  to  the  air  is  not  rate  limiting.  If  that 
had  been  the  case,  a  build-up  of  contaminant  would  have  been  expected  near  the 
surface.  This  was  not  observed  in  the  soil  gas  profiles  (Marrin  and  Thompson, 

1984: 173).  Therefore,  the  flux  given  in  the  study  for  the  California  site  will  be 
considered  the  flux  which  actually  transfers  to  the  air  from  the  soil. 

Transport  of  TCE  Through  the  Cron  Canopy.  The  development  below  seeks 
to  analytically  determine  the  expected  TCE  concentrations  from  the  assumed  level  of 
soil  gas  flux.  The  approach  will  be  to  first  develop  a  method  of  finding  the  effective 
difiiisivity  between  the  ground  surface  and  any  elevation  within  the  crop  canopy. 
With  this  difiiisivity,  the  concentration  gradients  between  the  ground  and  any 
elevation  can  be  found  from  Pick's  Law.  The  concentration  at  the  upper  boundary  of 
the  crop  canopy  is  assumed  to  be  zero.  With  the  concentration  at  the  top  zero,  the 
concentration  at  the  ground  surface  can  be  found  from  the  difference  in  the 
concentrations  at  the  two  elevations.  Calculation  of  the  concentration  differences,  or 
gradients,  between  the  ground  surface  and  other  elevations  is  then  used  to  find  the 
concentration  at  a  given  time  and  elevation. 

Difilisivitv  in  Com  Canopy.  The  difiusivities  within  a  com  canopy 
were  evaluated  by  Parmele  (1972).  The  vapor  transfer  coefficient  was  calculated  at 
the  top  of  the  com  over  a  24  hour  period  by  observation  and  calculation  of  the  flux 
of  water  vapor.  This  coefficient  was  then  distributed  with  decreasing  elevation 
through  the  crop  by  a  coefficient  which  expresses  the  difiiisivity  at  that  elevation  as  a 
fraction  of  the  one  at  the  top.  Because  the  v^or  transfer  coefficient  varies  with  time 
and  the  coefficient  which  distributes  it  through  the  various  elevations  varies  with 
hei^t,  the  diSusivity  varies  both  with  time  of  day  and  with  elevation  in  the  crop 
canopy. 
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The  distributing  coefiFiciait  will  be  referred  to  as  The  values  given  by 
Parmele  (1972)  are  shown  graphically  in  Figure  3  and  the  specific  data  points  are 
listed  in  ^pendix  1 . 


Figure  3:  Relative  dijffiisivity,  ip  (cnfi/s),  plotted  as  a  function  of  height,  x 

(Parmele  and  others,  1972:446) 

Though  Parmele  does  not  explain  the  zero  value  of  the  distributing  coefficient 
at  the  surface,  it  will  be  assumed  to  be  some  anomaly  of  the  calculations  he 
performed.  If  it  were  truly  representative,  it  would,  by  Pick's  Law,  represent  an 
in^jermeable  layer  at  the  surface.  At  the  surface,  the  diffiisivity  will  be  assumed  to 
equal  the  gaseous  molecular  diffusion  coefficient  for  TCE.  The  gaseous  diffusion 
coefficient  for  TCE  is  7030  cm^/day  (Peterson  and  others,  1988:576).  Values 
between  the  data  points  given  by  Parmele  will  be  calculated  by  linear  inteipolation 
using  the  linear  interpolation  function  of  Mathcad  edition  4.0.  These  interpolated 
values  will  be  represented  as  v|;(x). 

The  values  of  vapor  transfer  coefficient,  K2,  have  been  adapted  from  Parmele 
(1972:444).  All  of  the  variables  and  symbols  used  in  the  model  and  its  development 
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are  summarized  in  Appendix  4.  They  are  shown  graphically  in  Figure  4  and  the 
specific  data  points  given  by  Parmele  (1972)  are  listed  in  Appendix  2. 


Figure  4:  Vapor  transfer  coefficient,  (cm^/sec),  as  a  fiinction  of  time,  t 

(24  hr  clock) 

The  vqwr  transfer  coefficient  is  assumed  to  be  zero  during  most  of  the  night 
hours.  The  diffiisivity  at  these  times  virill  be  assumed  in  the  model  to  be  equal  to  the 
gaseous  diffiision  coefficient  for  TCE,  as  this  would  seem  to  represent  the  situation 
reasonably  even  under  the  most  stable  atmospheric  conditions.  The  stable 
conditions  of  the  night  atmosphere  represent  a  time  when  molecular  diffiision 
dominates  the  physical  transport  mechanisms. 

The  diffiisivity  at  a  specific  elevation  and  time  can  be  found  by  multiplying 
the  vapor  transfer  coefficient  for  time  of  day  by  the  distribution  coefficient  for  that 
elevation.  The  vapor  diffiisivity  at  the  top  of  the  canopy  at  time  t  will  be  represented 
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as  K^O).  These  values  will  be  interpolated  between  data  points.  KA,  the  difliisivity 
at  a  specific  elevation  and  time,  can  be  evaluated  as  follows; 

KA=  K.(t)H>(x) 


Effective  Difilisivitv.  A  point  at  any  elevation  sees  a  net  diffusivity 
between  itself  and  the  ground  surface.  Regardless  of  the  fluctuations  in  diflusivity 
with  elevation,  the  cc.icentration  at  the  point  is  only  a  functional  of  the  effective 
difiusivitv(KE)  between  that  elevation  and  the  ground.  In  the  problem  being 
considered,  the  diffusivity  varies  continuously  with  elevation.  This  continuous 
variation  is  viewed  here  as  a  series  of  constant  diffusivities  in  incremental  layers.  By 
analogy  to  hydraulic  conductivity  in  saturated  groundwater  flow,  the  effective,  or 
equivalent,  diffusivity  of  these  layers  can  by  written  as  follows  ; 


KB(x)  =  - 

y 


KE(x)=the  effective  diffusivity  at  elevation  x 
X  =  the  elevation  being  evaluated 
Axj  =the  incremental  layer  assumed  at 
increment  i 

KAi  =  the  diffusivity  which  is  assumed 

constant  through  the  incremental  layer, 
Axj  (Ritzi,  1993) 


Liss  (1974: 181)  has  confirmed  the  validity  of  the  analogy  to  groundwater  flow  and 
the  use  of  the  equation. 

The  development  of  the  effective  diffusivity  assumes  that  the  system  is  at 
steady  state.  With  this  assumption,  the  mass  flux  into  the  system  is  necessarily  the 
same  as  the  mass  flux  out  of  the  system  and  the  mass  flux  at  any  level. 
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Evaluation  of  TCE  Concentrations  in  the  Crop  Canopy.  The  diffusion  of  the 
soil  gas  through  the  crop  canopy  is  a  unique  situation.  The  source  of  TCE  is 
considered  essentially  infinite.  For  modeling  the  above  ground  e3q)osure,  the  soil  gas 
will  be  assumed  to  be  a  constant  source. 

Pick's  Law  can,  assumining  steady  state  difffusion,  be  written  (Couliette, 

1994): 

F  =  -KA^ 

Ax 


AC  =  the  change  in  concentration  over  Ax 
Ax  =  the  distance  over  which  the  change 
in  concentration  is  being  evaluated 


Rearrangement  of  Pick's  Law  gives  the  following; 


AC=- 


FAx 

KA 


With  the  assumption  that  the  concentration  is  zero  at  the  top  of  the  com,  the 
change  in  concentration  from  the  top  of  the  com  to  the  bottom  of  the  com  can  be 


calculated.  This  is  written  as  follows: 


Csgo  = 


FX 

KT 


Csgo=the  TCE  concentration  at  ground 
surface 

KT  =  the  effective  diffusivity  between  the 
top  of  the  com  and  the  ground 
X  =the  height  of  the  com,  assumed  to  be 

240cm 


Effective  diffiisivities  are  evaluated  at  various  heights  in  the  com  and,  based 
on  those  difhisivities,  the  steady  state  concentration  profile  at  a  given  point  in  time 
can  be  calculated.  The  concentration  at  elevation  within  the  canopy  can  be  written, 
by  once  again  rearranging  Pick's  Law,  as  follows; 
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Csg.  =  Csgo  + 


Fx 

KE. 


Csgx=the  TCE  concentration  at  xf 
X  =  the  elevation  being  evaluated 

KEx  =the  effective  diffusivity  for  level  x^ 


By  substituting  into  the  above  equation  for  Csgx ,  the  model  of  equilibrium 
concentration  profiles  due  to  soil  gas  at  time  t  can  be  written: 


FX 


Ft 


—  +  - 

KT 

X 

*  Ay 

L^ke(I)J 

Irrigation  Activity 

Calculation  of  Water  Requirements.  The  calculation  of  crop  water 
requirements  is  based  on  a  paper  by  Doorenbos  and  Pruitt  (1977).  All  the  necessary 
data  has  been  obtained  for  Hill  AFB  and  is  summarized  in  Table  3.  The  month  of 
July  is  assumed  for  the  time  of  year. 

The  initial  step  in  the  calculation  process  offered  by  Doorenbos  and  Pruitt  is 
the  estimation  of  evapotranspiration.  Of  the  methods  discussed  in  the  paper,  the 
Blaney-Criddle  was  selected  because  the  meteorological  parameters  it  uses  are 
readily  available  and  the  circumstances  at  Hill  AFB,  such  as  latitude  and  geography, 
meet  the  criteria  established  in  the  text.  The  relationship  is  as  follows: 

ETo=a[p(0.46T  +  8)] 

ETo  =  evapotranspiration  coefficient  in  mm/day 
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T  =  average  temperature  in  ®C 
p  =  average  percentage  of  yearly  daytime  hours  which 
occur  during  the  time  under  consideration,  tabulated 
by  Doorenbos  and  Pruitt  based  on  latitude  and  time  of 
year 

G  =  empirical  adjustment  factor  based  on  climate 
(D^enbos  and  Pruitt,  1977:3) 


Parameter 

Symbol 

Value 

Source 

Average  Temperature 

T 

24.4  ®C 

Intemat  Station  Met 
Climate  Summ 

Hill  AFB  Latitude 

41  North 

Intemat  Station  Met 
Climate  Summ 

Average  Relative 
Humidity 

RH 

28.4 

Intemat  Station  Met 
Climate  Summ 

Average  Rainfall 

Pe 

1.4  cm 

Intemat  Station  Met 
Climate  Summ 

Mean  Wind  Speed 

7.6  KTS 

Intemat  Station  Met 
Climate  Summ 

Table  3:  Hill  AFB  climatological  data. 


The  evapotranspiration  coefficient  incorporates  the  influence  of  climate.  All  climate 
data  was  acquired  from  the  CD  ROM  International  Station  Meteorological  Climate 
Summary  which  gives  data  from  the  weather  station  at  Hill  AFB.  Due  to  the 
empirical  adjustment  factor,  g,  Doorenbos  and  Pruitt  (1977)  separate  out  the  value, 
f,  as  follows: 

f  =  [p(0.46T+8)] 

The  evapotranspiration  coefficient  is  then  found  based  on  data  tabulated  by 
Doorenbos  and  Pruitt  (1977)  based  on  minimum  relative  humidity,  daylig^  hours, 
and  wind  estimates.  The  final  number  will  be  obtained  from  tabulated  data 
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(Doorenbos  and  Pruitt,  1977:fig  1)  assuming  medium  RH,  mean  daytime  wind  of  5 
m/s,  and  medium  to  high  sunshine  hours  (Doorenbos  and  Pruitt,  1977:3). 

Another  coefficient,  the  crop  coefficient  incorporates  the  influence  of  the 
specific  crop.  The  crop  coefficient  used  here  assumes  sweet  com  as  would  more 
likely  be  grown  in  a  private  garden.  Together  the  two  coefficients  are  used  to 
calculate  the  crop  evapotranspiration.  The  relationship  is  as  follows: 

ETcrop  =  KcETo 


Kc  =  crop  coefficient 
ETcrop  =  crop  evapotranspiration 


The  net  irrigation  requirement  is  given  by  the  following  relationship: 

In  =  30-  ETcrop-  (Pe+  Ge+  Wb) 


In 

net  irrigation  requirement  (mm/month) 

Pe  = 

average  rainfall  (mm/month) 

Ge  = 

groundwater  contribution  (mm/month) 

Wb  = 

water  stored  in  the  soil  at  the  start  of  the  time 
period  (mm/month) 

In  the  model,  the  water  stored  in  the  soil  and  the  groundwater  contribution  will  be 
assumed  to  be  zero.  The  total  volume  required  for  irrigation  is: 


'lO' 

In 

.Ep, 

a-LR)_ 

V  =  total  volume  of  irrigation  required  (mm/month) 

Ep  =  percent  efficiency  of  the  irrigation  method 

LR  =  leaching  requirement  (fraction) 

This  equation  is  slightly  different  than  the  one  given  by  Doorenbos  and  Pruitt 


(1977).  Doorenbos  and  Pmitt  multiply  the  net  irrigation  requirement  by  the  acreage 
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to  get  a  cubic  volume.  Leaving  this  out  gives  millimeters  of  water,  which  is  better 
suited  to  this  model. 

The  irrigation  system  will  be  assumed  to  be  100%  efficient.  The  leaching 
requirement  coefficient  increases  the  irrigation  volume  so  that  salts  are  leached 
away.  The  parameters  necessary  to  calculate  LR,  such  as  "electrical  conductivity  of 
the  soi  .ration  extract  for  a  given  crop  appropriate  to  the  tolerable  degree  of  yield 
reduction,"  (Doorenbos  and  Pruitt,  1977:77)  are  not  available.  The  tabulated  data 
in  the  text  has  been  used  to  arrive  at  an  assumed  value  of  0. 1 49.  This  is  an 
approximate  upper  bound  on  LR  and  therefore  maximizes  V.  Maximizing  V 
maximizes  the  mass  of  TCE  applied  to  the  crop,  and,  therefore,  is  a  conservative 
estimate. 


Parameters 

Symbol 

Value 

Source 

Percentage  Yearly 
daytime  hours 

P 

33% 

(Dorrenbos, 
1977:Table  1) 

Crop  Coefficient 

Kc 

1.05 

(Dorrenbos, 
1977:Table  21) 

Table  4;  Parameters  in  model  of  irrigation  requirements. 


Losses  During  Spray  Irrigation.  One  method  of  delivery  of  water  modeled 
was  spray  irrigation.  The  model  treated  the  aeration  of  the  water  during  irrigation  as 
a  major  loss  of  TCE  to  the  troposphere.  Any  TCE  volatized  before  reaching  the 
ground  was  assumed  to  be  unavailable  to  the  com  crop.  Giardino  (1992)  examined 
three  models  of  losses  jfrom  a  shower  spray.  The  models  estimate  a  mass  transfer 
coefficient  and  use  that  coefficient  to  calculate  the  concentration  remaining  in 
solution.  The  model  used  here  was  chosen  because  more  of  the  parameters  utilized 
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are  available  and  because  in  the  experimental  study  by  Giardino  (1992)  the  error  it 
demonstrated  consistently  underestimated  volatilization  during  spraying.  This  will 
allow  the  model  to  remain  conservative.  The  other  models  examined  by  Giardino 
(1992)  occasionally  overestimated  the  volatilization. 

The  mass  transfer  coefficient  is  given  as  follows  (Giardino  and  others, 
1992:1602): 


,  .0037SV^ 

l+fp^/p.) 

=  the  liquid  phase  mass  transfer  coefficient 
Vj  =  the  velocity  of  the  water  drop 
mi=  the  viscosity  of  water 
=the  viscosity  of  air 


The  constants  used  in  the  calculation  are  summarized  in  Table  4.  Utilizing 


this  the  concentration 


remaining  in  solution  can  be  written  as  follows: 


C^i  =Coexpj 


Cli  =  the  concentration  of  TCE  after  it  reaches  the 
ground 

Cq  =  the  initial  concentration  of  TCE  in  the  irrigation 
water 

d  =  the  drop  diameter 

td  =  the  residence  time  of  the  drop  in  air 

This  is  slightly  different  from  the  equation  given  by  Giardino  (1992).  In  this 
model,  the  air  flow  will  be  assumed  to  be  sufficient  to  make  the  concentration  of 
TCE  surrounding  the  drops  negligible.  This  assumption  is  justified  by  a  shower 
study  by  McKone  and  Knezovich  (1991)  where  the  transfer  efficiency  did  not 
change  significantly  with  shower  length  (McKone  and  Knezovich,  1991 :835). 
Longer  showers  would  have  had  increasing  gaseous  concentrations  of  TCE  around 
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the  drops.  Based  on  this  assumption,  the  gaseous  concentration  surrounding  the 
drops  in  the  original  model  given  by  Giardino  was  set  to  zero  and  the  equation 
written  accordingly. 

Some  of  the  parameters  used  in  the  spray  model  will  have  to  be  assumed. 

The  study  presenting  the  model  is  for  a  shower  spray,  and  the  drop  diameter  and 
residence  time  will  differ  from  that  of  spray  irrigation  equipment.  The  drop  diameter 
will  be  assumed  to  be  .25  cm.  This  is  within  the  range  for  which  the  model  is  still 
valid  but  larger  than  most  of  the  drop  sizes  observed  by  Giardino  (1992).  The 
longer  trajectory  required  of  a  spray  irrigation  system  would  be  consistent  with  a 
larger  drop  size(Giardino  and  others,  1992: 1603).  Larger  drops  also  tend  to  be  less 
efQcient  at  volatilizing  the  TCE,  and  it  is  therefore  a  conservative  assumption 
(Gillham  and  Rao,  1990:493).  The  velocity  is  assumed  to  be  507  cm  /  s.  This  is  the 
velocity  found  by  Giardino  (1991 )  for  a  10  L  /  min  shower  spray. 

The  residence  time  will  be  approximated  by  the  following  trajectory  equation 
(Serway,  1986:61): 

^  sin  (a) 

— 

g 

a  =  initial  angle  of  trajectory  (assume  30°) 
g  =  acceleration  due  to  gravity  (980  cm  /  s^) 

This  is  a  rough  estimation  where  the  frictional  forces  due  to  air  have  been 
neglected.  Based  on  the  values  already  assumed,  t(}  is  calculated  to  be  .26  sec. 

Losses  Prior  to  Flood  Irrigation.  The  irrigation  water  would  likely  be 
transported  to  the  com  crop  by  flow  through  a  culvert.  Any  TCE  which  volatilizes 
wfrile  in  the  culvert  is  assumed  to  be  lost  to  the  troposphere  where  it  cannot  expose 
the  com  crop.  One  method  of  estimating  the  rate  of  volatilization  is  to  estimate 
exponential  decay  rate  of  the  TCE  concentration  in  the  aquifer.  Based  on  that  rate, 
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Parameter 

Symbol 

Value 

Source 

Viscosity  of  Water 

1.0  X  10"3  N-sec  / 
m2 

(Serway, 

1986:333) 

Viscosity  of  Air 

Ha 

1.78  X  10^  N-sec  / 
m2  (at  15  “C) 

(Kuethe  and 
Chow,  1986:518) 

Table  5:  Constants  used  in  spray  model. 


the  amount  of  time  that  the  water  takes  to  reach  the  field  can  be  used  to  directly 
calculate  the  amount  of  TCE  which  remains  when  the  water  reaches  the  crop. 


Chiou  (1979)  gives  a  direct  formulation  of  the  mass  transfer  rate.  The  mass  transfer 


coefficient  is  given  by: 


k  =  op  K 


,  2  uRT  j 


k  =  mass  transfer  coefficient  from  water 
(cm/s ) 

Kh=  Henry’s  law  constant  (dyne-cm/g) 

M  =  molecular  weight  of  the  solute  (g/mole) 

T  =  temperature  (®K) 

P  =  coefficient  which  is  a  function  of  air  flow 
(unitless) 

a  =  coefficient  which  is  a  function  of  the 
liquid  turbulence  (unitless) 

R  =  universal  gas  constant  (J  /  mole  °K) 


Chiou  evaluates  a  and  p  for  a  number  of  situations.  Converting  the  units  of 
the  Henry's  constant  for  TCE  gives  1 .07  x  1 0^  dyne-cm/g  under  the  following 
assumptions: 


-Kh  for  TCE  is  1 . 1 8  kPa  m^/mol  at  25®C 


-dyne  =  1  g  cm  /  s2  (Weast,  1985;F-77) 

-1  Pascal  =  1  N  /  m2  (Serway,  1986:317) 

For  this  Henry's  constant  and  unstirred  water,  Chiou  indicates  an  a  of  .9  (Chiou, 
1979:235).  The  a  for  unstirred  water  is  used  under  the  assumption  that  the  flow  in 
the  culvert  is  laminar  or  non-turbulent.  Using  the  a  for  turbulent  flow  would  give  a 
larger  evaporation  rate  and  would  therefore  be  a  less  conservative  assumption. 

Chiou  (1979)  also  indicates  a  P  for  TCE  evaporating  into  still  air  as 
1.83  X  10~^.  The  p  for  still  air  is  used  so  that  the  estimate  is  conservative  under 
even  the  most  stable  atmospheric  conditions.  Both  of  these  coefficients  are  unitless. 
The  transfer  coefficient  can  be  used  as  a  decay  rate  and  the  calculation  of  a 
concentration  is  as  follows  (Masters,  1991:72): 


Co=  the  initial  concentration  of  TCE 
t  =  the  time  that  the  water  takes  to  reach  the 
field 

Time  in  the  equation  can  be  replaced  by  the  distance  traveled  divided  by  the 
velocity  of  the  water.  The  equation  used  in  the  model  will  be: 

f  -U  ^ 

1  =  the  distance  traveled  from  well  to  field 

Vw=  the  velocity  of  the  water  in  the  culvert 

The  velocity  of  the  water  will  arbitrarily  assumed  to  be  1  m/s,  and  the  distance  will 

be  assumed  to  be  1000  meters. 
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rate  dependent  on  the  initial  soil  conditions,  and  decreases  to  the  rate  of  the  saturated 
hydraulic  conductivity.  This  initial  infiltration  rate  is  not  currently  well  defined  and 
a  review  of  the  literature  has  not  found  an  analytical  model  to  portray  it.  This  model 
will  assume  that  the  infiltration  rate  is  constant  at  the  saturated  hydraulic 
conductivity.  The  saturated  hydraulic  conductivity  will  be  assumed  to  be  lO’^ 
(cm/sec).  This  is  the  low  end  of  the  hydraulic  conductivities  that  might  be  expected 
for  a  sihy  sand  such  as  are  present  at  Hill  AFB.  According  to  Darcy's  Law,  the 
average  linear  velocity  under  these  conditions  will  be: 


In  this  case,  hydraulic  potential,  dh/dz,  is  the  depth  of  the  ponded  water.  The 
variables  K  and  n^  are  the  hydraulic  conductivity  and  the  soil  porosity  respectively. 
The  soil  porosity  will  be  assumed  to  be  .6,  which  is  a  typical  number  for  a  sihy  sand 
(Fetter,  1988:67). 

In  the  model  for  the  penetration  of  the  contaminant,  the  amount  of  time 
during  which  the  contaminated  irrigation  water  is  infiltrating  the  surface  represents 
the  parameter  tQ.  The  model  will  assume  weekly  irrigation.  This  can  be  evaluated 
as  follows: 

Vper 


The  variable  Vper  is  the  depth  of  water  applied  during  the  weekly  irrigation  and  is 
sin:^)ly  one  fourth  of  the  monthly  net  requirement.  The  rate  of  water  ^plication  is 
assumed  to  be  much  greater  than  the  rate  of  infiltration. 


Analytical  Estimation  of  Penetration  and  Distribution.  Models  describing  the 


leaching  of  organic  chemicals  were  designed  for  use  with  agricultural  chemicals. 
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They  assume  a  finite  amount  of  surface  applied  chemical  being  leached  by  a  large 
quantity  of  uncontaminated  water.  The  case  of  contaminated  irrigation  water  is 
unique  in  that  the  amount  of  organic  applied  is  proportional  to  the  amount  irrigation 
water.  The  model  used  here  assumes  that  a  chemical  dissolved  in  water  enters  the 
soil  from  time  zero  to  to  and  that  after  to  pure  water  enters  the  soil.  In  the 
development  of  the  model,  the  time  to  represents  the  time  when  the  chemical  on  the 
surface  is  depleted. 

The  model  was  first  developed  by  Oddson  (1970).  The  development  and  re- 
evaluation  of  the  equations  is  given  by  Letey  and  Oddson  (1972).  The  depth  into  the 
ground  is  z.  The  coordinate  system  used  has  ground  surface  as  z=0  and  increasing  z 
into  the  ground. 

Letey  and  Oddson  (1972)  present  two  separate  equations.  The  adsorbed 
concentration  as  a  function  is  given  by ' 

S(z,t)=K^Ci,iexp 
The  boundary  conditions  are: 

c(0,t)  =  CLi  forO:St<to 
c(0,t)=0  fort>to 

These  conditions  state  that  a  solution  of  concentration  Cl]  is  applied  until  time  tQ 
and,  after  that  time,  the  concentration  applied  is  0.  Iq  is  the  modified  Bessel 
function. 

A  condition  for  this  specific  equation  is  that; 

S(z;t)  =  0  for  0  <  t  <  z  /  Vz 

This  condition  simply  states  that  the  adsorbed  concentration  at  point  z  is  zero  before 
the  wetting  front  reaches  it. 


-oK. 


f) 

j  exp(-T)Io 


oK^zt 


th 
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The  result  of  this  equation  gives  the  adsorbed  concentration,  in  the  same  units 
of  Cli ,  at  any  given  time  and  point.  The  equation  is  not  valid  for  time  after  the 
infiltration  of  contaminated  water  has  stopped.  Letey  and  Oddson  give  a  modified 
form  of  the  equation  for  time  after  to  which  will  be  used  to  investigate  the  impact  of 
rain.  The  development  of  the  equation  assumes  that  pure  water  is  entering  the  soil 
after  to  as  would  be  the  case  for  rain. 

The  coefficients  a  and  are  the  adsorption  rate  constant  and  the  adsorption 
equilibrium  constant  respectively.  is  readily  available  and  is  listed  Table  1  for 
TCE.  The  transformation  of  to  the  K<}  used  in  the  equation  is  direct  and  will  be 
discussed  below.  Care  must  be  taken  as  it  is  dependent  on  some  of  the 
characteristics  of  the  soil.  The  value  used  here  is  one  experimentally  determined  by 
Paviostathis  and  Jagial(  1991).  It  has  the  advantage  that  the  soil  it  is  based  on  is 
well  described  in  the  study  and  the  soil  is  reasonably  close  to  that  of  Hill  AFB. 
Several  other  soil  parameters  will  be  taken  from  the  study  by  Paviostathis  and  Jagial 
(1991).  Taking  the  parameters  from  an  ^:tual  soil  has  the  advantage  of  knowing 
that  the  parameters  assumed  do  coexist  in  an  actual  soil  and  therefore  do  not  conflict. 

The  a  describes  the  rate  at  which  the  equilibrium  between  the  soihed  and 
solution  concentrations  reach  equilibrium.  Mathematically  it  can  be  stated  as: 


—  =  a(K,C-S) 

at  ' 


S  is  the  sorbed  concentration  and  K  is  the  partition  coefficient  at  equilibrium.  At 
equilibrium  the  rate  of  change  of  the  sorbed  concentration  goes  to  zero  as  S  =  KpC. 
Currently  sufficient  data  is  not  available  to  calculate  a  specific  a  for  TCE  in  a  given 
soil.  None  of  the  sorption  studies  found  evaluate  the  necessary  rate  information. 
Letey  and  Oddson  did  provide  a  range  of  values,  .033  to  5.09,  which  have  been 
estimated  for  various  soils  and  organic  chemicals.  This  range  will  be  used  to 
observe  the  influence  on  the  depth  and  concentration  of  the  penetrating  contaminant. 
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It  is  necessary  to  clarify  some  of  the  other  parameters  used  by  Letey  and 
Oddson  (1972)  and  Oddson  (1970)  as  they  differ  slightly  from  the  traditional  ones; 

Vz  -  The  velocity  Vi  given  by  the  Darcy's  Law  and  saturated 
hydraulic  conductivity  is  normally  used.  In  this  case  Vz  is  the 
velocity  in  the  pore  space  given  by: 

Vp=Vz/e 

6  is  the  volumetric  water  content  and  will  be  assumed  to  be  22%, 
which  is  the  value  established  for  the  soil  being  used  by 
Paviostathis  and  Jagial  (1991 :276). 

S  -  The  quantity  of  adsorbed  material  is  normally  given  on  a 
mass/mass  basis.  The  relationship  between  the  traditional 
quantity.  S',  and  the  S  used  in  the  model  is  as  follows: 

s  =  ps' 

P  is  the  bulk  density.  The  bulk  density  of  the  soil  in  the  study  by 
Pavlotathis  (1991)  is  calculated  to  be  2.8  g  /  cm^.  Thus,  the 
product  this  equation  is  a  mass  of  contaminant  per  cm^  of  soil. 

Kp  -  The  normal  partitioning  coefficient,  K^,  is  related  to  the  one 
utilized  in  the  model  by; 

Kp=  KdP/e 

Cg  -  The  normal  solution  concentration,  C,  is  given  in  mass  per 
cm^  of  solution.  The  C  used  in  the  model  is  related  by  the 
following: 
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Cs-ec 

The  result  is  a  mass  of  contaminant  in  liquid  phase  per  volume  of 
soil. 


In  the  model,  all  of  the  values  are  calculated  in  the  traditional  units. 

The  liquid  phase  concentration  at  lime  after  to  is  as  follows  (Letey  and 
Oddson,  1972:411): 

8(2, t) 


C.(z,t)  = 


K. 


+  Cl,  exp 
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This  m.odel  does  make  several  assumptions  which  keep  it  from  being  truly 
predictive.  One  problem  is  the  assumptions  it  makes  about  adsorption  and 
desorption.  The  coefficient  a  is  only  valid  if  the  rate  of  sorption  is  constant  over 
time.  This  may  not  be  true  for  TCE,  and  it  is  definitely  not  true  for  some  organic 
chemicals.  The  model  also  assumes  that  sorption  is  entirely  reversible,  and  this  may 
also  be  a  faulty  assumption  (Oddson  and  others,  1970:416).  In  addition,  it  neglects 
the  vjqxir  phase  interactions  which  may  be  present.  The  current  state  of  research  is 
not  able  to  describe  the  adsorption  behavior  of  vapor  phase  TCE  in  soils  (Peterson 
and  others,  1988:571).  The  error  involved  in  not  including  the  vapor  phase  may  or 
may  not  be  significant  because  at  the  time  the  penetration  is  being  modeled  most  of 
its  pore  space  will  be  occupied  by  solution.  The  model  is  still  a  useful  description  of 
chemical  behavior  and  provides  insights  into  the  factors  that  influence  it. 
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Exposure  In  the  Air  Phase  Due  to  Irrigation 

For  modeling  the  air  under  the  crop  canopy,  the  water  reaching  the  crop  will 
be  assumed  to  be  deposited  immediately  under  the  crop  and  pond  there.  For  the 
evaluation  of  the  exposure  in  the  air  under  the  com  canopy,  the  water  will  be 
assumed  to  not  infiltrate  the  soil.  This  will  tend  to  overestimate  the  exposure  that 
would  actually  take  place  as  this  makes  more  mass  ofTCE  readily  available.  In 
reality,  the  when  the  TCE  infiltrates  the  soil  it  would  evaporate  much  more  slowly 
due  to  sorption  and  limited  exposure  to  the  air. 

The  model  used  will  treat  the  ponded  TCE  aqueous  solution  as  a  homogenous 
concentration  which  interfaces  with  the  air  at  x  =  0.  The  coordinate  system  will 
assume  that  x  is  positive  in  the  vertical  up  direction.  All  of  the  water  will  be 
assumed  to  be  deposited  at  time  t  =  0. 

The  difiusivity  varies  both  with  elevation  and  with  time  of  day.  Solutions  to 
the  difiusion  differential  equations  which  allow  this  dual  variation  have  not  been 
located.  However,  a  solution  which  allows  the  difiusivity  to  vary  over  time  is 
available.  The  concentration  will  be  evaluated  at  individual  elevations  using  the 
effective  difiusivity  discussed  earlier.  With  the  elevation  held  constant,  the 
difiusivity  varies  only  with  time.  These  points  can  be  assembled  to  obtain  a 
complete  concentration  profile  for  the  com  canopy.  Crank  (1970:12)  presents  a 
solution  for  a  planar  source.  The  solution  is  as  follows: 

C(x,,k)  =  — ,  exp 

The  evaluated  X  is  substituted  to  allow  the  difiusivity  to  vary  over  time.  Crank 
(1970)  gives  the  variable  as  follows: 

t 

X  =  jK^(x)ih 
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One  unique  aspect  of  this  physical  situation  is  the  air/water  interface.  There 
is  a  transfer  coefficient  in  the  liquid  and  gas  film  layer  which  must  be  included  in  the 
evaluation  of  Kg. 
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IV.  Results 


Air  Exposure  Due  to  Soil  Gas. 

The  crop  canopy  concentrations  (ppm)  of  TCE  resulting  from  a  flux  of 
12  p,g/m2-day  are  plotted  for  three  different  times  of  day  in  Figures  6  through  8. 
The  quantity  of  a  i^g/cm^  has  been  converted  to  ppm  on  a  volume/volume  basis 
assuming  that  the  unit  of  1  ^g^cm^  is  equal  to  2  ppm  ( ATSDR,  1 989:33).  The 
difiusivity  values  which  lead  to  these  concentrations  are  shown  in  Figure  S  for  0800 
and  1200.  The  molecular  difiusion  that  is  dominant  during  the  night  hours  is 
negligibly  small  in  comparison  to  these  values  and  was  thus  not  shown.  All  of  the 
graphs  in  this  section  are  summarized  in  Table  6. 


Figure  5:  Difiusivity  (cm/s)  as  a  function  of  elevation  (cm)  within  the  crop 
canopy  at  0800  and  1200  hours. 
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Figure  6:  Concentration  (ppm)  of  TCE  at  elevation  (cm)  within  the  crop 
canopy  at  0200  hours  due  to  a  flux  of  12  ^g/m2-day. 


Concentration 

(ppm) 


Elevation 

Icm) 


Figure  7:  Concentration  (ppm)  of  TCE  at  elevation  (cm)  within  the  crop 
canopy  at  0800  hours  due  to  a  flux  of  1 2  p-g/m^-day. 
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Figure 

Distance 

Flux  Rate 
(pg/m2-day) 

TimeofDav 

(hours) 

Parameter 

Modeled 

5 

Elev.  =  0-240 
(cm) 

12 

0800, 1200 

Dififtisivity 

(cm/s) 

6 

Elev.  =  0-240 
(cm) 

12 

0200 

Concentration 

(ppm) 

7 

Elev.  =  0-240 
(cm) 

12 

0800 

Concentration 

(ppm) 

8 

Elev.  =  0-240 
(cm) 

12 

1200 

Concentration 

(ppm) 

Table  6:  Summary  of  graphs  used  to  describe  soil  gas  activity. 
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Losses  During  Transport. 

Spray  Irrigation.  The  spray  of  a  TCE  contaminated  water  is  modeled  in 
Figure  9  below.  The  concentration  arriving  at  the  crop  is  shown  as  a  function  of  the 
concentration  in  the  original  aquifer  water.  The  TCE  which  is  lost  to  the  atmosphere 
is  considered  to  be  removed  and  has  no  potential  to  expose  the  com  crop.  The 
concentrations  are  given  in  p.g/L  and  therefore  correspond  to  ppb  on  a  mass/volume 
basis. 

0.08 

0.06 

Concentration 
at  the  Crop 

0.(H 

0.02 

0 

Figure  9 .  The  concentration  (ng/L)  of  TCE  reaching  the  crop  due  to  spray 
irrigation  is  shown  as  a  function  of  the  initial  concentration  of  TCE  in  the 
aquifer  (pg/L). 

Flood  Irrigation.  The  concentration  of  TCE  in  the  water  arriving  at  the  com 
crop  during  flood  irrigation  is  modeled  in  Figure  10  below.  The  amount  of  TCE 
delivered  is  shown  as  it  varies  with  the  distance  that  the  water  has  to  travel  to  get  to 
the  crop.  The  model  assumes  that  the  TCE  is  transported  in  a  culvert  with  a  uniform 
depth  of  .5  m  and  that  the  initial  concentration  is  1 6  p-g/L. 
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15 
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at  Crop 

Vaialion  With  DiitaMc 


0 

Figure  10:  The  concentration  (jiig/L)  of  TCE  reaching  the  crop  due  to  flood 
irrigation  is  shown  as  a  function  of  the  distance  (m)  the  TCE  travels  to  the 
field  in  transport. 

Crop  Exposure  Due  to  Irrigation 

Air  Exposure.  The  following  figures  show  the  results  for  the  transport  of  the 
TCE  up  through  fh©  crop  canopy.  Figures  1 1  throu^  13  model  the  variation  of 
concentration  with  elevation  and  time  after  a  0600  irrigation.  As  discussed 
previously,  the  units  of  air  concentration  are  ppm  on  a  volume/volume  basis.. 

The  concentration  is  next,  in  Figure  14,  modeled  at  an  elevation  of  10  cm  for 
an  extended  period  after  irrigation  to  observe  the  attenuation  of  the  TCE.  In  Figure 
1 5,  the  attenuation  is  observed  during  the  extreme  outlying  hours.  All  of  the  graphs 
in  this  section  are  summarized  in  Table  7. 

Exposure  In  Soil.  In  Figure  16,  the  aqueous  concentration  of  TCE  is  shown 
as  a  ftinction  of  the  depth  at  the  time  when  the  water  completes  the  infiltration  into 
the  soil.  The  units  of  liquid  phase  TCE  are  ng'L,  which  corresponds  to  roughly  a 
ppb  on  a  mass/mass  basis.  In  Figure  1 7,  the  concentrations  of  sorbed  and  liquid 
phase  TCE  are  shown  in  units  of  ^g/cm^.  Assuming  a  bulk  density  of  2.8  g/cm^. 
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1  fig/cm^  corresponds  to  approximately  .36  ppb  on  a  mass/mass  basis. 


Concentration 

1  Bout 


6  Hoiin 


Figure  1 1 :  The  concentration  of  TCE  (ppm)  is  plotted  as  a  function  of  the 
elevation  within  the  crop  canopy  (cm)  at  a  time  1  hour  and  6  hours  after  a 
0600  irrigation.  The  data  assumes  16  pg/L  concentration  reaching  the  crop 
and  1 .4  cm  of  average  rainfall. 


Figure  12:  The  concentration  of  TCE  (ppm)  is  plotted  as  a  function  of  the 
elevation  within  the  crop  canopy  (cm)  at  a  time  1  hour  and  12  hours  after  a 
0600  irrigation.  The  data  assumes  16  pgT  concentration  reaching  the  crop 
and  1.4  cm  of  average  rainfall. 
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Figure  13:  The  concentration  of  TCE  (ppm)  is  plotted  as  a  function  of  the 
elevation  within  the  crop  canopy  (cm)  at  a  time  1  hour  and  30  hours  after  a 
0600  irrigation.  The  data  assumes  16  ^g/L  concentration  reaching  the  crop 
and  1.4  cm  of  average  rainfall. 


Time 


Figure  14:  The  concentration  of  TCE  (ppm)  at  a  height  of  10  cm  is  plotted 
as  a  function  of  time  (hours)  after  a  0600  irrigation.  The  data  assumes 
16  |xg/L  TCE  reaching  the  crop  and  1.4  cm  of  average  rainfall. 
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Figure  15:  The  concentration  of  TCE  (ppm)  at  a  height  of  10  cm  is  plotted 
as  a  function  of  time  (hours)  after  a  0600  irrigation.  The  data  assumes 
16  \x%fL  concentration  reaching  the  crop  and  1 .4  cm  of  average  rainfall. 


Figure 

Time 

After 

Irrigation 

(Hours) 

Concentration  in 
Irrigation  Water 
(lig/L) 

Elevation 

(cm) 

Parameter 

Modeled 

11 

1,6 

16 

0-240 

12 

1,12 

16 

0-240 

Concentration 

(ppm) 

13 

1,30 

16 

0-240 

14 

0-500 

16 

10 

15 

400-700 

16 

10 

IBSSH 

Table  7:  Summary  of  graphs  used  to  model  the  crop  air  exposure. 
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The  aqueous  concentration  of  TCE  is  modeled  as  it  varies  with  the 
parameters  of  sorption  rate  and  sorption  equilibrium.  In  Figure  18,  is  held 
constant  and  the  impact  of  a  variation  of  a  is  observed.  In  Figure  19,  a  is  held 
constant  and  the  concentrations  are  observed  for  two  different  values  of  K^.  All  of 
the  graphs  used  to  model  the  exposure  in  the  soil  phase  are  summarized  in  Table  8 
below. 

Affect  of  Rainfall.  The  influence  of  rainfall  is  modeled  as  pure  water  being 
applied  after  the  contaminated  water  has  completely  entered  the  soil.  Figure  20 
illustrates  the  influence  of  the  rainfall  on  the  concentration  profile  of  the  liquid  phase 
contaminant. 


Figure  16:  The  concentration  of  TCE  (jig/L)  in  aqueous  phase  in  the  soil  is 
plotted  as  a  function  of  depth  (mm).  The  data  assumes  t  =  100  sec, 
concentration  reaching  the  crop  is  16  p-g/L,  a  =  .05,  and  =  1 1.7  ml/g. 
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Depth 

Figure  17:  The  mass  of  TCE  ()xg/cm^)  sorbed  (solid  line)  and  in  the  aqueous 
phase  (pg/cm^)  (dashed  line)  is  plotted  as  a  function  of  depth  (mm).  The 
data  assumes  t  =  100  sec,  K<j  =  1 1.7  ml/g,  a  =  .05,  and  the  concentration  of 
TCE  reaching  the  crop  is  16  pg/L. 


Figure  18:  The  concentration  of  TCE  in  aqueous  phase  (pg  /  L)  at  a  =.0 1 
(solid  line)  and  the  concentration  of  aqueous  TCE  (pg  /  L)  at  a  =.  1  (dashed 
line)  are  plotted  with  depth  (mm).  The  data  assumes  t  =100  sec  and  K^= 
11.7  ml/g. 
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Depth 


Figure  19:  The  concentration  of  TCE  (fig  /  L)  in  aqueous  phase  at  Kj  = 
11.7  ml/g  (solid  line)  and  at  K<}  =  6.35  ml/g  (dashed  line)  are  plotted  as  a 
function  of  depth  (mm).  Data  assumes  t  =  100  sec,  a  =  .05  ml/g,  and 
concentration  reaching  the  crop  is  16  )ig  /  L. 


Aqueous  Concentration 


Figure  20:  The  concentration  of  TCE  in  the  aqueous  phase  (|ig/L)  is  shown 
as  a  fimction  of  depth.  The  data  assumes  t  =  200  sec.  Kd=  11.7  ml/g,  a  =  .05 
and  the  concentration  of  TCE  reaching  the  crop  is  16  iig/I  . 
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Figure  Sorption  Equil. 
Coeflf. 
(ml/g) 


Concentration 

Time 

in  Water 

After 

(lig^) 

Irrigation 

(sec) 

16 

100 

1 

16 

100 

16 

100 

16 

100 

16 

1 _ 1 

200 

Aqueous 

Cone. 

(^g/L) 


Aqueous 

and 

Sorbed 

Cone. 


Aqueous 

Cone. 


Aqueous 

Cone. 

(^g/L) 


Aqueous 

Cone. 

(Mg/L) 


Table  8:  Summary  of  graphs  used  in  modeling  the  leaching  of  ^ '  *  'nto  the  soil 
environment. 
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V.  Discussion  and  Conclusions 

Air  E?cposure  Due  to  Soil  Gas 

The  air  phase  exposure  due  to  soil  gas  is  shown  in  Figures  6, 7  and  8.  At 
night  the  concentrations  are  expected  to  be  highest  due  to  lower  diffusivity.  The 
concentrations  are  all  parts  per  trillion  in  magnitude.  The  concentrations  during  the 
more  turbulent  daytime  hours  are  all  more  than  4  orders  of  magnitude  less  than  those 
at  nig^t  and  would  quite  likely  be  undetectable.  It  is  important  to  remember  that  this 
model  was  done  with  the  assumption  that  the  system  had  reached  steady  state.  The 
actual  numbers  are  possibly  quite  lower  than  the  ones  calculated  here,  but  it  is  safe 
to  say  that  they  are  not  higher.  The  levels  of  TCE  in  the  parts  per  trillion  range, 
lacking  significant  biomagnification  in  the  com,  are  considered  insignificant  and  the 
actual  numbers  are  likely  much  lower  than  those  calculated.  The  exposure  due  to 
soil  gas  is  therefore  considered  negligible  in  the  air  phase. 

Transport 

Spray  Irrigation.  The  modeling  of  transport  of  water  to  the  crop  by  spray 
irrigation  is  shown  in  Figure  9.  Based  on  the  model  used,  more  than  99%  of  the 
TCE  can  be  expected  to  be  removed  from  the  water.  While  this  varies  significantly 
from  the  study  conducted  by  Giardino  (1992),  the  study  by  Giardino  did  not  assume 
that  the  ambient  concentration  was  negligible  as  was  done  in  this  model.  The  study 
by  Giardino  (1992)  was  examining  volatilization  in  a  shower  and  the  assumption 
was  therefore  not  appropriate.  The  99%  removal  of  TCE  provides  an  opportunity  to 
reduce  the  exposure  of  the  com  by  selection  of  irrigation  systems. 
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Flood  Irrigation.  The  losses  during  flood  irrigation  are  shown  in  Figure  10. 
The  volatilization  using  this  method  of  irrigation  is  much  lower  and,  for  the  numbers 
that  would  be  expected  reasonable  in  an  agricultural  situation,  no  significant  loss  of 
TCE  occurs.  For  the  losses  to  be  significant,  the  distance  traveled  by  the  water  to 
the  crop  has  to  be  on  the  order  of  hundreds  of  kilometers. 

Crop  Exposure  Due  to  Irrigation 

Air  Exposure.  The  exposure  in  the  air  phase  is  shown  in  Figures  1 1  through 
13.  The  concentrations  are  all  less  than  10  ppm.  These  concentrations  are  very  near 
concentrations  which  can  be  breathed  directly  without  significant  risk  to  human 
health  (ATSDR,  1989:3).  Without  significant  biomagnification,  the  exposure  can 
be  considered  negligible.  It  is  noteworthy  that  these  concentrations  are  only  present 
in  the  very  lowest  regions  of  the  crop  canopy.  Figure  14  shows  a  profile  over  time 
for  a  height  of  10  cm.  This  Figure  demonstrates  that,  even  as  low  as  10  cm,  the 
concentrations  do  not  exceed  a  ppm.  The  height  of  the  concentration  is  significant 
because  leaves,  where  transpiration  tak^  place,  are  the  most  likely  place  for  uptake 
of  TCE  (Bleckmann,  1994).  The  higher  concentrations  are  much  less  likely  to  be 
available  for  uptake.  In  addition,  the  assumption  of  ponding  is  conservative  and  the 
leaching  of  contaminated  water  would  limit  the  contamination  to  levels  below  the 
ones  found  here. 

Soil  Exposure.  The  aqueous  TCE  concentration  is  shown  in  Figure  16.  The 
level  of  contamination  below  2  mm  is  shown  to  be  in  the  ppb  range.  The  majority  of 
the  mass  of  TCE  is  retained  in  the  top  2  mm  of  soil.  These  results  are  consistent 
with  those  of  Letey  and  Oddson  (1972).  At  these  levels  of  soil  penetration  it  is 
doubtful  that  the  roots  would  be  exposed.  Figure  1 7  compares  the  sorbed  mass  to 
the  aqueous  mass.  The  importance  of  this  is  that,  as  mentioned  previously,  the  sorbed 
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mass  is  not  available  for  uptake  by  the  plant.  The  model  demonstrates  that  nearly  all 
of  the  mass  of  the  TCE  is  sorbed  by  the  soil  and  therefore  not  going  to  be  taken  up 
by  the  plant. 

The  impact  of  the  sorption  parameters  is  modeled  in  Figures  18  and  19.  If 
either  the  rate  of  sorption  or  the  sorption  equilibrium  coefficient  decreases 
significantly  the  depth  of  penetration  by  the  aqueous  TCE  increases.  If  it  is  not 
sorbed  to  the  soil  rapidly  or  if  the  soil  capacity  to  soi  _  l.  is  low,  the  TCE 
continues  to  flow  into  the  lower  levels  of  the  soil  where  the  roots  may  be  exposed.  In 
a  moist  soil  environment,  the  sorption  is  mostly  by  way  of  the  organic  matter  in  the 
soil.  The  organic  matter  in  the  soil,  therefore,  plays  an  important  role  in  determining 
IQj  and,  hence,  the  depth  of  penetration  by  the  TCE.  Very  little  is  known  about  the 
sorption  rate  coefficient  for  TCE.  In  order  to  get  a  good  idea  of  the  actual 
penetration  more  research  will  have  to  be  conducted  to  better  define  this  parameter 
for  TCE.  From  the  results  of  the  model,  it  can  be  concluded  that  the  properties  of  the 
soil  are  important  in  determining  the  likely  exposure  to  TCE. 

Figure  20  shows  the  impact  of  rain  on  the  TCE.  The  TCE  is  pushed  down 
through  the  soil  in  a  wave,  and  the  more  rain  takes  place  the  more  TCE  reaches  the 
area  of  the  roots.  The  rain  therefore  increases  the  likelihood  of  the  uptake  of  the 
TCE  by  the  plant. 

The  rainfall  amounts  modeled  in  Figure  20  were  approximately  twice  the 
average  monthly  rainfall  at  Hill  AFB.  Based  on  the  model  in  Figure  24,  a  com  crop 
outside  Hill  AFB  could  reasonably  exp^  to  have  200  ppb  or  less  of  aqueous  TCE 
at  a  depth  of  less  than  4  mm.  While  200  ppb  may  be  a  potentially  risky  amount  of 
TCE  if  it  is  available  for  uptake  by  the  plant,  the  depth  of  penetration  is  still  minimal 
and  will  mitigate  the  effects.  It  is  important  to  note  the  interesting  result  that 
irrigation  with  16  ppb  aqueous  solution  of  TCE  resulted  in  a  contamination  at  the 
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200  ppb  level.  This  may  be  due  to  the  sorption  properties  of  the  soil  retarding  the 
penetration  of  the  TCE  relative  to  the  penetration  of  the  water. 

Conclusions 

Based  on  a  review  of  the  literature,  a  methodology,  consisting  primarily  of  a 
series  of  models,  was  developed  for  evaluation  of  the  exposure  pathways  for  TCE 
between  a  contaminated  aquifer  and  a  food  crop.  Using  information  gathered  about 
the  climate  and  environmental  conditions  near  Hill  AFB,  the  model  was  used  to 
examine  the  potential  for  a  problem  at  that  location  due  to  an  existing  TCE 
contamination.  The  methodology  did  provide  significant  information  about  the 
potential  for  a  risk  to  human  health.  The  most  important  contributions  are  as 
follows: 

-the  air  exposure  due  to  soil  gas  is  negligible 
-the  method  of  irrigation  is  an  important  factor  in  the  exposure 
-the  air  exposure  due  to  irrigation  is  negligible  and  much  of 
the  resulting  vapor  phase  TCE  is  unavailable  for  uptake  by 
the  plant 

-the  soil  exposure  is  mitigated  by  the  lack  of  penetration  of  the 
TCE  and  the  potential  for  exposure  is  influenced  by  the 
properties  of  die  soil 

-rainfall  assists  the  penetration  of  the  TCE  toward  the  roots 

The  methodology  does,  therefore,  contribute  useful  information  about  factors 
influencing  risk  to  human  health.  With  the  types  of  information  that  are  provided, 
the  risk  analysis  and  prevention  at  a  contaminated  site  can  be  guided  much  more 
effectively. 
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Appendix  1 :  Distribution  coefficient  data  points  (Parmele  and  others,  1 972:446) 


Elevation  in  Cron  Canoov,  x 
(cm) 

Distribution  Coefficient  u/ 
(unitless) 

240 

1.0 

227 

.85 

189 

.50 

135 

.25 

81 

.17 

37 

.02 

0 

_ 

0 

68 


Appendix  2:  Vapor  transfer  coefficients  at  the  top  of  the  com  (Parmele  and  others, 

1972:444) 
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Appendix  3:  Parameters  and  variables  used  in  the  literature  review 


Variable 

Svmbol 

[  Units 

Octanol/Water  Partition 
Coefficient 

^Q)W 

Unitless 

Octanol/Carbon  Partition 
Coefficient 

Koc 

Unitless 

Vapor  Pressure 

None 

1  ! 

mmHg 

Aqueous  Solubility 

None 

g/L 

Henry's  Law  Coefficient 

Kh 

Pa  m^/mol  or 

Adsorption  Equilibrium 
Constant  for  TCE 

Kd 

mL/g 

Molecular  Weight  TCE 

M 

g/mol 

Concentration 

C 

Diffiision  Coefficient 

D 

cm^/sec 

Flux 

F,Jv 

Porosity 

•I) 

Unitless 

Volume  of  Soil 

V 

cm^ 

Volume  of  Voids 

Vy 

cm3 

Tortuosity  Factor 

Sorbed  Mass  of 
Contaminant 

S 

lAg 

Absorption  Rate 

a 

Unitless 
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Appendix  4:  Parameters  and  variables  used  in  the  model 


Parameter  or  Variable  i 


Symbol 


Units 


Difllisivity  Distributing 
Coefficient 


Unitless 


Elevation  Within  Com 
Canopy 


Diffusivity  at  Top  of 
Com  Crop 


Diffiision  Coefficient  at 
Elevation  within  Canopy 


Effective  Diffusivity 


Effectivity  Diffusivity  at 
Top  of  Com 


Height  of  Com 


Concentration  of  TCE 
Reaching  Surface  Under 
Crop 


Soil  Gas  Flux 


Concentration  at  Ground 
Surface 


Average  Temperature 


Average  Rainfall 


Crop  Evapotranspiration 

ETo 

mm  /  day 

Leaching  Requirement 

LR 

unitless 

Efficiency  of  Irrigation 
System 

EP 

% 

Ground  Water 
Contribution 

Ge 

mm/day 

Water  Stored  In  Soil 

Wb 

mm/day 

Empirical  Correction 
Factor 

cr 

unitless 

mean  percentage  of 
annual  daytime  hours 

P 

% 

Uncorrected 

Evapotranspiration 

coefficient 

f 

mm/day 

Crop  Evapotranspiration 

ETcrop 

mm  /  day 

Crop  Evapotranspiration 
Coefficient 

Kc 

unitless 

Net  Irrigation 
Requirement 

In 

mm  /  month 

Total  Volume  of 
Irrigation  Water  Applied 

Vi 

mm  /  month 

Weekly  Irrigation 
Volume 


Vper 


mm  /  week 


Appendix  4  continued 


Liquid  Phase  Mass 
Transfer  Coefficient 

kl 

cm  /  s 

Velocity  of  Water  Drop 

Vd 

cm  /  s 

Viscosity  of  Water 

Nsec  /  nfi 

Viscosity  of  Air 

N'sec  /  rrfi 

Concentration  of  TCE  in 
Aquifer  Water 

Co 

Drop  Diameter 

d 

cm 

Residence  Time  of  Drop 
in  Air 

sec 

Initial  Trajectory  Angle 
of  Spray 

a 

degrees 

Mass  Transfer 
Coefficient  From  Water 

k 

cm/  sec 

Molecular  Weight  of 
TCE 

M 

g  /  mole 

Air  Flow  Coefficient 

P 

unitless 

Liquid  Tuihulence 
Coefficient 

a 

unitless 

Universal  Gas  Constant 

R 

J/mole  °K 

Velocity  of  the  Water  in 
the  Culvert 

Vw 

cm/s 
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Distance  From  Well  to 
Crop 

1 

meters 

Uniform  Depth  of 
Culvert 

L 

meters 

Conductivity 

K 

cm  /  sec 

Porosity 

ne 

unitless 

Initial  Depth  of  Water 
Applied 

Vi 

cm 

Time  for  Complete 
Infiltration 

to 

sec 

Velocity  of  Infiltrating 
Water 

Vz 

cm /sec 

Depth  into  the  Ground 

z 

cm 

Sorbed  Mass 

S 

pg/  cm^ 

Sorbed  Mass 

S' 

lig/g 

Concentration  of  TCE  in 
Solution 

Cs 

^g/L 

Soil  Volumetric  Water 
Content 

e 

% 

Soil  Bulk  Density 

P 

pg  /  cm^ 
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Appendix  5:  Mathcad  4.0  code  used  to  represent  soil  gas  exposure 


The  foflowing  enters  the  soil  gas  levels  from  the  analogous  California  site  into  a  vector  so  that  the 
rTOt  exposure  due  to  soil  gas  can  be  shown 


i  =0,1.4 


Depth.  =.305zsgj 


The  concentrations  are  in  ppb.  The  depths  are  in  feet 
and  converted  to  meters. 


The  folowing  Mathcad  4.0  template  models  the  concentration  levels  within  the  crop  canopy  at  any 
given  elevation  and  time  due  to  a  given  soil  gas  flux. 


Hi(x)  -linterpfx.v.x) 


The  actual  data  points  from  Parmele  (1972)  and  the 
corresponding  elevation  within  the  com  canopy  are  entered 
Into  two  vectors. 


A  linear  interpolation  is  used  to  generate  a  distribution 
coefficient  which  is  a  continuous  function  of  height. 


The  water  vapor  diffusivity  at  the  top  of  the  com  as  adapted 
from  Pannele  ( 1972)  is  entered  into  a  vector  with  another 
vector  corresponding  to  the  time  of  day  for  the  data  point. 
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K(t)  =bnteip(£.K,t) 


Linear  interpolation  is  used  to  generate  diffusivites  which 
are  a  continuous  function  of  time. 


F  =  .0012  The  assigned  fhjx  is  entered  in  ug/day  per  square 

centimeter. 

X  =  1.4  ..  240  The  vanable  x  gives  the  elevation  in  1cm  increments,  up  to 

240  cm. 

A  X  ~  1 

The  increments  of  x  are  1  cm  in  width. 


KA(xt)  =iH:K(t)\»)(x)>0.K(t)H<x).. 001356) 


KE(x,t)  := - 

X 

s 

i=  1 


X 


Ax 

KA(i.l) 


K(x.t)  generates  a  diffusivity  in  cm2/sec  If  the 
ctiffusivity  from  Parmele  is  zero,  the  molecular 
diffusion  coefficient  for  TCE  is  retimed.  The 
equivalent  diffusivity  is  then  calculated  by 
evaluating  the  contributions  of  the  individuai 
layers. 


C«g0(t)  = 


F.J..J-.J-.240 
24  60  60 


K£(240.t) 


The  concentration  at  the  surface  is  calculated  in  ug/tm3. 


C8g(xt)  =Csg0(t)- 


- X 

24  60  60 

KE(xt) 


The  concentration  at  any  given  elevation  at  any  given  time  is 
calcLteted. 


The  following  portion  of  the  templet  evaluates  the  concentration  at  a  given  elevation  over  a  24hr 
period  or  due  to  a  range  of  soil  gas  fluxes. 


F  =0. 10..  100 
l  =0,0100..  2400 


The  elevation,  x,  must  be  specified  in  cm.  The  range  of  flux 
must  be  specified  in  ug/day  per  square  cm. 


CmFOIF,  t)  := —  The  concentration  at  the  surface  is  calculated  in  ug/cm3. 

KE(240.t) 

The  folowing  portion  of  the  templet  evaluates  the  concentration  at  a  given  elevation  over  a  24hr 
period  or  due  to  a  range  of  soil  gas  fluxes. 


F  =0, 10.100 
t  =0,0100..  2400 


C8gF0(F.t) 


F-J-  -L  J-240 
24  60  60 


KE(240,t) 


The  elevation,  x.  must  be  specified  in  cm.  The  range  of  flux 
must  be  specified  m  ug/day  per  square  cm. 


The  concentration  at  the  surface  is  calculated  in  ug/cm3. 


C8gF(F.t.x)  :=C8gF0(F.t) - ^  ^ 

KE(x,t) 


The  coTKentration  at  any  given  elevation  at  any 
given  time  is  calculated. 
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Appendix  6:  Mathcad  4.0  code  used  to  model  transport  and  irrigation  volume 


The  following  models  the  concentration  of  TCE  reaching  the  crop  from  spray  irrigation. 

g;:980  a=30  Vd  =507  d:=.23  The  assumed  values  for  the  initial  trajectory,  a.  the 

acceleration  due  to  gravity.g  (cm/s2).  the  drop 
CO  -  0. 10..  10000  d  =  0, 1000..  1000000  diameter,  d,  and  the  assumed  velocity. 


Id  = 


g 


td  =0.259 


The  residence  time  of  the  drops  in  air  is  calculated. 


pw  =1.0  10"^ 


pa  =  1.78-10* 


The  viscosities  of  water  and  air  are  given  in  kg^s. 


KL  = 


.00375Vd 


1  + 


pw 

pa 


The  mass  transfer  coefficient  for  the  drops  Is 
calculated. 


CSpi(C0)  :=C0  exp 


-KL 


The  coTKentration  of  TCE  in  the  irrigation  water  when 
it  reaches  the  crop  is  calculated  in  the  units  of  the 
aquifer  contamination. 


The  following  models  the  concentration  of  TCE  reaching  the  crop  due  to  flood  irrigation. 
M  =130  a  =.9 

P:=1.8310'’  KH  =1.0710* 


T  =24,4 


R  =8.13 


Kexch  :=a  P  KH' 


M 


1000 


2  ji  R100*  (T -I- 273.16) 


The  exchange  constant  is  evaluated. 


Vw.-i  L -.5  The  velocity  of  water  flow  In  the  culvert  (mfe)  and  depth  of  the 

culvert  (m)  are  given. 


CFl(C0,d)  :=  co-exp 


Kexch- d  \ 
.5-lOO-Vw/ 


The  exchange  constant  is  used  to  model  an  exponential 
decay  of  the  concentration  of  TCE  in  the  water. 


The  following  calculates  the  water  required  per  weekly  irrigation  as  a  function  of  the 
average  rainfall. 

p  :=  .33  T  := 24.4  Kc  - 1.05  The  average  dally  fraction  of  daytime  hours,  leaching 

requirement,  crop  coefficient,  and  average  temperature  in  July 
(degrees  C)  are  0ven. 


The  variable  f  is  used  to  find  the  crop  evapotranspiration,  ETo, 
from  the  Tabulated  data  (Doorenbos,  1977;fig  f ) 
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From  the  tabuteited  data  Ge -o  Wb  =o  Ep  - 1 

ETo  =8 

ETcrop  =  Kc  ETo  The  crop  water  requirement  is  calculated 


In(Pe)  =  30  ETcrop  -  (Pe  +  Ge+  Wb ) 


The  net  irrigation  reqiirement  is  calculated  as  a  function 
of  the  average  rainfall  The  groundwater  contribution.  Ge, 
and  the  water  initially  stored  in  the  soil  have  been  set  to 


Pe  -0. 1..  100 


zero. 


V(Pe)  =-L.i2^ 
Ep  1-  LR 


The  monthly  volume  of  water  rerpred  is  calculated.  The 
irrigation  system  efficiency.  EP.  is  assumed  to  be  100%  and 
the  leaching  requirement  has  been  set  to 


Vp«(Pe) 

4 


The  depth  in  mm  of  a  weekty  irrigation  event  is  calculated. 


78 


Appendix  7:  Mathcad  4.0  code  modeling  the  soil  and  air  exposure  due  to  irrigation 

The  folowing  Mathcad  4.0  template  models  the  exposure  due  to  spray  irrigation. 


0 

0 

■  0 

0  ' 

37 

,02 

0 

02 

81 

.17 

0 

04 

I  = 

135 

V  - 

.25 

300 

06 

189 

.5 

750 

08 

227 

.85 

1500 

10 

240 

1.0 

K  = 

1500 

£  = 

12 

«(<X)  = 

linterp(x.»ii.x) 

1500 

14 

The  actual  data  points  from 

700 

16 

Parmele  (1972)  and  the 

500 

18 

corresponrtng  elevation  within 

250 

20 

the  com  canopy  are  entered 

into  two  vectors. 

0 

22 

A  linear  interpofeitlon  is  used  to  generate  a  distribution 

0 

24. 

coefTicierrt  which  is  a  continuous  function  of  height 
The  waiter  vapor  diffusivity  at  the  top  of  the  corn  as  adapted 
from  Parmele  (1972)  is  entered  into  a  vector  with  another 
vector  conesponding  to  the  time  of  day  for  the  data  point. 


K(t)  :=linletp(e,K.t) 


Linear  interpolation  is  used  to  generate  diffusivities  which 
are  a  continuous  function  of  time. 


t;=l,10.40 


Time  is  entered  in  hours  for  use  in  exponential  decay. 


The  variable  x  gives  the  elevation  in  1  cm  increments,  up  to 
240  cm. 

The  increments  of  x  are  1  cm  in  width. 

The  following  calculates  the  water  required  per  weekly  irrigation  as  a  function  of  the 
average  rainfall. 

p  =  .33  T  .=  24.4  Kc  =1.05  The  average  daily  fraction  of  daytime  hours,  leaching 

requirement,  crop  coefficient,  and  average  temperature  in  July 
(degrees  C)  are  given. 


x:=l,4..240 

iix:=l 


f  =p(.46T  +  8) 

From  the  tabulated  data. 
ETo  =8 

ETcrop  .=  Kc  ETo 


The  va;  lable  f  is  used  to  find  the  crop  evapotranspiration,  ET o, 
from  the  Tabidated  data  (Doorenbos.  1977;fig  1 ) 

Ge  =0  Wb  -0  ‘ 

The  crop  water  reqiarement  is  calculated. 


79 


Iii(Pe)  =30  ETcrop  -  (Pe  +  Ge-t-  Wb  ) 


The  net  irrigation  requirement  is  calculated  as  a  function 
of  the  average  rainfall  The  groundwater  contribution,  Ge, 
and  the  water  initaHy  stored  in  the  soil  have  been  set  to 


zero. 


V(Pe) 

Ep  1  -  LR 

Vpe<Pe) 

4 


The  monthly  volume  of  water  required  is  calculated  The 
irrigation  system  efficiency.  EP,  is  assumed  to  be  100%  and 
the  leaching  requirement  has  been  set  to 

The  depth  in  mm  of  a  weekly  irrigation  event  is  calculated 


The  following  models  the  diffusion  into  the  crop  canopy  of  the  TCE  which  is  ponded  at 
the  base  of  the  corn. 


M  130  o  :=  .9  The  values  for  molecular  weight  (g/mol),  alpha,  beta,  Henry's 

-1S3  10'^  KH  -1 07  10*  constant  (dyne-cm/g),  temperature(degrees  K).  and  ideal 

gas  constant  (J/mole-K)  are  given 

T  :=24.4  R  :=8.13 


Kint  Is  the  exchange  constant  for  the  water/air 
interface  and  is  of  the  form  required  to  account  for 
the  contribution  of  the  interface  to  the  net 
diffusivity. 


KA(x.t) 


ii:  K(mo<l(  t.  24)  )i(<x)>0.  K(  mod(  t.  24)  )i|)(x),  ,001356) ,  iffKf  t)- 1*<  X)  >0,  K(t  )•»(<  X) .  .001356) 


KE(x.t)  := - - 

(_L),  y 

\Kint/  ^  KA(i,t) 
i=  1 


KA(x,t)  generates  a  cSffusivity  in  cm2feec.  If 
the  diffusivity  from  Parmele  is  zero,  the 
molecular  diffusion  coefficient  for  TCE  is 
returned.  The  equivalent  diffusivity  is  then 
calculated  by  evaluating  the  contributions  of  the 
individual  layers. 
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F(CO,t)  -Kiiu  lcO—^  expl 
\  1000  ' 


I  Kint  t-60  60\\ 

\  .5100  II 


Di(x.t) 


ft 


J6 


KE  k - di 

3600/ 


The  variable  of  substitution  and  the  diffusion 
equation  con^lete  the  modeling 


CLl  =o..5..  10 

11  =0, 10. 500 

12  ^  400.410..  700 

MASS(CL2.Pe)  =  VpeifPe)  10  CL2 .001 


The  mass,  MASS,  of  TCE  deposited  on  the  ground 
per  unit  area  is  determined  This  calculation 
assumes  that  the  initial  aqueous  concentration 
reaching  the  com  Is  given  in  ug/L  (ppb).  The  volume 
of  irrigation  is  converted  from  mm  to  cm  and  ugA.  is 
converted  from  L  to  cubic  centimeters 


The  concentration  in  ug/cm3  is  evaluated 
as  function  of  time  (Crank, 1970:10). 

X  MASS(CL2.Pe) 

CA(CL2,Pe.xt)  =  ■  ^  '  e: 

^■Di(x3600(t+6)) 

The  following  tiKxlels  the  leaching  of  the  chemicat  into  the  ground. 


4Di(x3600-(t  +  6)) 


e:=.22  con  =10'^  Pe  =140  B:=2.8  z:=.5,1..7 

CL2  =1,10.100  a::.05..1..  1  Kd  =1,3.20 


Vi  :=  c<mVpa(Pe)- 10 


^  _Vp«(Pe) 

Vi 

Vi  B 

V2:=—  Kp(Kd):  =  -Kd 

6  6 


The  velocity  of  the  infiltration  is  calculated  based  on  the  volume 
of  irrigation  and  conckjctvity  of  the  soil. 

The  time,  tO,  for  all  of  the  water  to  infiltrate  the  soil  is  calculated 
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at - 

Vz 


S(z,l,a.Kd.CL2)  -Kp(Kd)  CL2  6  .OOl  exp  a  Kp{Kd) — 

\  Vz, 


exp(  t)  Id  2 


a  Kp(Kd)  z  T 


SS(z,t.a.Kd.CL2)  =  -  S(z.l.a.Kd,CL2) 
B 


The  concentration  is  formulated  in  ug  TCE  per  gram 
of  soil 


The  amount  of  TCE  in  the  liquid  phase  is  calculated. 


S(z,t.a.Kd.CL2) 

C(z,t.tt.Kd.CL2)  - -  +  i 

Kp(Kd) 


/-aKp(Kd)z\ 

CL26-.001exp - — — —  c: 

\  Vz 


4H- 


n  ‘-- 

— ]  10  2  Kp{Kd)  z  -^ — — 

Vz/  Vz 


CS(z,t,a.Kd,CL2)  :=-^z,t.a,Kd.CL2)  1000  The  concentration  is  formulated  in  ug  TCE  per  Liter 
^  of  solution. 

The  folowing  models  the  effect  of  pure  water  being  applied  to  the  crop  after  irrigalion. 

The  following  models  the  effect  of  pure  water  being  appNed  to  the  crop  after  irrigation. 

Craii<z.t,a.Kd.CL2)  ^ -n 

Kp<Kd) 


t-— 1  t-— -100 

10  2-  a*  Kp(Kd)  z  -^ — ^  -  exp(a  l00)  I0  2  a*  Kp(Kd)  z  -^ — — - ^ 

.  L  I  Vz  J  J  [  [  Vz 


CR(z.t,a,Kd,CL2)  :=-^H(z,t,a.Kd.CL2) 

e 
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